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ABSTRACT 

We have analysed radial velocity measurements for known transiting exoplanets to 
study the empirical signature of tidal orbital evolution for close-in planets. Compared 
to standard eccentricity determination, our approach is modified to focus on the re- 
jection of the null hypothesis of a circular orbit. We are using a MCMC analysis of 
radial velocity measurements and photometric constraints, including a component of 
correlated noise, as well as Bayesian model selection to check if the data justifies the 
additional complexity of an eccentric orbit. We find that among planets with non-zero 
eccentricity values quoted in the literature, there is no evidence for an eccentricity 
detection for the 7 planets CoRoT-5b, WASP-5b, WASP-6b, WASP-lOb, WASP-12b, 
WASP-17b, and WASP-18b. In contrast, we confirm the eccentricity of HAT-P-16b, 
e = 0.034 ± 0.003, the smallest eccentricity that is reliably measured so far for an 
exoplanet as well as that of WASP-14b, which is the planet at the shortest period 
(P = 2.24 d), with a confirmed eccentricity, e = 0.088 ± 0.003. As part of the study, 
we present new radial velocity data using the HARPS spectrograph for CoRoT-1, 
CoRoT-3, WASP-2, WASP-4, WASP-5 and WASP-7 as weU as the SOPHIE spectro- 
graph for HAT-P-4, HAT-P-7, TrES-2 and XO-2. 

We show that the dissipative effect of tides raised in the planet by the star and 
vice- versa explain all the eccentricity and spin-orbit alignment measurements available 



for transiti ng planets. We revisit the mass-period relation (Mazeh et al. 2005 Pont 
et al.|[20TT | and consider its relation to the stopping mec hanism of orbital m igration 



for hot Jupiters. In addition to CoRoT-2 and HD 189733 (iPont et al.|2009b|), we find 



evidence for excess rotation of the star in the systems CoRo i'-18, HA'i'-P-20, 
and WASP-43. 
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1 INTRODUCTION 

Most of the information we have about the formation, evo- 
lution and structure of exoplanets have come from the study 
of transiting planets. This is possible because the combina- 
tion of radial velocity measurements with transit photom- 
etry can provide powerful constraints on the physical and 



* Based on observations made at the 1.93-m telescopes at Ob- 
servatoire de Haute-Provence (CNRS), Prance with the SOPHIE 
spectrograph. 



orbital parameters of an exoplanet, such as the planetary 
mass, radius, orbital eccentricity, etc. 

A selection effect due to geometry means that most 
transiting planets with radial velocity confirmation are 
found on very short period orbits with P ~ 1 — 20 days. 
The close-in planets with periods of a few days are expected 
to experience strong tidal effects (e.g. |Rasio fc Ford||1996[ ), 
which should increase sharply with decreasing period and 
these orbits are thus expected to circularise on a timescale 
much smaller than the system age. A higher tendency for 
such circular orbits is indeed observed in the sample of tran- 
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siting planets, as compared to those from radial velocity 
surveys. This has been interpreted as a signature for tidal 
circularization. The transition from eccentric orbits to cir- 
cular orbits at short period has also been seen in binary star 



0.5). 



(1) 



systems, e.g. Mathieu & Mazeh ( 1988 \ and Mazeh ( 2008 1 



Over the last few years, we have carried out a monitoring 
programme to obtain several radial velocity measurements 
of known transiting planetary systems with the intention of 
refining the orbital properties such as orbital eccentricity 
and spin-orbit alignment angle. We have used the SOPHIE 
spectrograph in the Northern hemisphere and the HARPS 



spectrograph in the Southern hemisphere (e.g. Loeillet et al, 
2008; H ebrard et al. 2008: Husnoo et al.||2011| |Pont et"aL 
2011, ESO Prog. 0812.C-0312). 



One issue is the difficulty of measuring the orbital eccen- 
tricity of exoplanets for faint stars, especially for low-mass 
planets. While it is impossible to prove that an orbit is cir- 
cular, with e = exactly, we can place an upper limit on the 
eccentricity of a given orbit (e.g. reject e > 0.1 at the 95% 
confidence level). In fact, a number of eccentric orbits have 
been detected at short period, but follow-up observations us- 
ing photometry or additional radial velocity measurements 
led to the conclusion that some of these eccentricities had 
originally been overestimated. For example, the WASP- 10 
syste m (Christian et al. 2009} was revisited by Maciejewski 
et al. (2011b), who showed that the initially reported eccen- 
tricity (e = O.O59IQ QJ4) had been overestimated and was in 
fact compatible with zero. The orbital eccentricity of WASP- 
12b | |Hebb et al.||2009 



e = 0.049 ± 0.015) is similarly com- 



patible with zero ( Husnoo et al.|[2011| ), and the original de- 
tection was possibly due to systematic effects (weather con- 
ditions, instrumental drifts, stellar spots or scattered sun- 
light). 

The eccentricity distribution at short period has a cru- 
cial importance for any theory of planetary formation and 
orbital evolution. Planets on orbits that are consistent with 
circular gather in a well-defined region of the mass-period 



plane, close to the minimum period for any given mass ( Pont 
|et al.|pC)ll| ). We now show that there are no exceptions to 
this pattern, and revisit some apparent exceptions as re- 
ported in the literature. As an ensemble, the totality of tran- 
siting planets considered in this study are in agreement with 
classical tide theory, with orbital circularisation due to tides 
raised on the planet by the star and tides on the star raised 
by the planet, to varying degree depending on the position 
of the planet-star system in the mass-period plane. 

In this study, we consider new radial velocity measure- 
ments made with HARPS and SOPHIE, as well as mea- 
surements present in the literature. We use photometric 
constraints in the form of the orbital period P and mid- 
transit time Ttr, both of which can be measured accurately 
using transit photometry, and we also consider constraints 
from the secondary eclipse where available. In fact, if we 
define the orbital phase <j) = (t — To)/P to be zero at mid- 
transit time To — Ttr, a planet on a circular orbit would 
have a mid-occultation phase of <^ = 0.5 (by symmetry). A 
planet that is on an eccentric orbit will, however, have a 
mid-occultation phase different from 0.5 (unless the orbital 
apsides are aligned along the line of sight). This allows us 
to place a constraint on the ecoscj projection of the eccen- 



to first order in e, where we now define (jtocc to be the 
phase difference between the mid-transit time and the mid- 
occultation time, i.e. 0occ = {Tsec ~ Ttr) / P , where T^^^ is 
the time of the secondary eclipse following the transit time 
Ttr- The component e sin lj is dependent on a ratio involv- 



ing the durations of the occultation and transit ( Winn et al. 
|2005| ), 



Ttra 



T 



Ttra 4" Toe 



(2) 



to first order in e, where Ttra and Tocc are the transit and 
occultation durations respectively, although this constraint 
is weaker than the one on ecosw. 

In addition to the reanalysis of radial velocity measure- 
ments with photometric constraints, we also introduce two 
modifications to the Markov Chain Monte Carlo (MCMC) 
process commonly used by teams analysing radial velocity 
data to work out the orbital parameters of transiting ex- 
oplanets. This involves a new treatment of the correlated 
noise present in most radial velocity datasets, as well as 
analysing the data in model selection mode to check if an 
eccentric orbit is indeed justified, given the additional com- 
plexity of the eccentric version of a Keplerian orbit. 

The present time is significant in the study of exoplan- 
ets, because a number of high quality measurements are now 
available for the three main observable effects of tides: circu- 
larisation, synchronisation and spin-orbit alignment. In Sec- 
tion |2] we describe our new radial velocity measurements 
obtained with SOPHIE and HARPS for 10 objects, as well 
as the measurements we collected from the literature for 
this study. We then describe the analysis we performed, 
in Section [S] In Section |4] we describe the objects in the 
classes "eccentric", "compatible with circular, e < 0.1", and 
"poorly constrained" and present the updated eccentricities, 
as shown in Table [Tg] In Section [5j we consider our orbital 
eccentricity estimates in the light of tidal effects inside the 
planet due to the star and vice-versa. We find that our re- 
sults are compatible with classical tidal theory, removing 
the need for perturbing stellar or planetary companions to 
excite non-negligible eccentricities in short period orbits. 



Winn et al. (20101 presented a discussion of the avail- 



able measurements of the projected spin-orbit alignment an- 
gles, and found that hot planet-hosting stars (Tgfj > 6250 K) 
had random obliquities whereas cooler stars (Tgg- < 6250 K) 
tended to have aligned rotations. These authors suggested 
that this dichotomy can be explained if all these stars har- 
bouring a planetary system start off with a random obliquity 
following some dynamical interaction, but only cool stars 
with a significant convective layer are able to undergo tidal 
effects leading to alignment. We verify in Section [5] that the 
strong exceptions WASP-8 and HD 80606 are indeed sys- 
tems with weak tidal interactions, and that the observation 
that these two are misaligned, is not incompatible with tidal 
theory. 



In a number of cases, such as HD 189733 (Henry & 
Wlmi|[2008t, WA SP- 19 ( |Hebb et al.][20Tol ) and CoRoT-2 



(Lanza et al. 2009 1, the rotational period of the star is known 



tricity, as given by Winn et al. ( 2005 I (slightly modified) 



from photometric monitoring. Assuming the results of |Winn| 
|et al.| ( |2010| | are correct in the sense that the convective 
layer in G dwarfs would cause tidal dissipation that aligns 
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the stellar equator with the planetary orbit, the negligible 
value of the projected spin-orbit angle A means that the 
obliquity is indeed zero, i.e. the stellar equators are aligned 
with the orbital planes. In this case, a measurement of the 
projected equatorial rotational velocity of the star (iisini) 
through Doppler broadening yields the rotational period of 
the star. This means that for G dwarfs at least, we have 
enough information to observe the effect of tidal interac- 
tions on the stellar rotation. We show in Section |5] that in 
addition to CoRoT-2 and HD 189733 ( |Pont et al.|2009b| ), we 
find evidence for excess rotation of the star in the systems 
CoRoT-18, HAT-P-20, WASP-19 and WASP-43. 

The preliminary results from this study were published 
in Pont et al. (20111, and the exact numerical values of 



the eccentricities have been updated in this paper to reflect 
our new choice of radial velocity measurement correlation 



timescale r — 1.5 d (see Section 3.21, as opposed to r = 0.1 
d in Pont et al. (2011 1. The overall results, i.e. the clear sep- 



aration between orbits that are consistent with circular and 
eccentric orbits in the mass-period plane, does not change in 



this paper. The mass-period relation of Mazeh et al. ( 2005 1 
is still clearly present, with low-mass hot Jupiters on or- 
bits that are consistent with circular clumping in a defi- 
nite region of the mass-period plane, with heavier objects 
moving closer in, to shorter periods. This strongly suggests 
that tidal effects are involved in the stopping mechanism 
of these objects. A similar stopping mechanism can be seen 
at higher planetary masses, but destruction of the planet 
is not excluded in many cases. Other effects, such as spin- 
orbit alignment and stellar spin-up also point strongly to- 
wards the scenario of 'Rasi o fc Ford| |1996) where the short 
period orbits of hot Jupiters are formed by dynamical scat- 
tering, which produces eccentric and misaligned orbits. This 
is followed by tidal dissipation which leads to circularisation 
at short period, spin-orbit alignment and synchronisation of 
the rotation of the host star. 



2 OBSERVATIONS 

We include 73 measurements for 6 objects with the HARPS 
spectrograph (Tables [T] [5] [s] [2] |5] and [6] ) and 45 measure- 
ments for 4 objects with the SOPHIE spectrograph (Ta- 



bles [7[[8j[9| and 10 1. Both are bench-mounted, fibre-fed spec- 
trograph built on the same design principles and their ther- 
mal environments are carefully controlled, to achieve pre- 
cise radial-velocity measurements. The two instruments have 
participated in the detection and characterisation of nu- 
merous transiting exoplanets, notably from the WASP and 
CoRoT transit searches. The wavelength calibrated high- 
resolution spectra from the instruments are analysed using 
a cross correlation technique which compares them with a 
mask consisting of theoretical positions and widths of the 
stellar absorption lines at zero velocity ( |Pepe et al.|2002| ). 

We carried out a literature survey and collected radial 
velocity measurements for 54 transiting planets, as well as 
other relevant data such as the orbital periods and the time 
of mid-transit. For the cases of CoRoT-1, CoRoT-2 and GJ- 
436, we also used the secondary eclipse constraint on the 



eccentricity component ecosw from Alonso et al. (20091, 



Alonso et al. (2009b I and Deming et al. (20071, respectively. 
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Table 1. HARPS radial velocity measurements for CoRoT-1 (er- 
rors include random component only). 
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4768.52895 


-56.5039 
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4769.51404 
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Table 2. HARPS radial velocity measurements for CoRoT-3 (er- 
rors include random component only). 



ing exoplanets, we had to stop the clock somewhere, and 
we picked the 1st of July 2010. We selected only objects 
that had been reported in peer-reviewed journals or on the 
online preprint archive ArXiV.org. Moreover, we selected 
systems with well measured parameters (planetary radius 
Rp and mass Mp to within 10%) and excluded faint objects 
iy > 15). At that time, 64 such systems were known. We 
reanalyse the existing radial velocity data for 54 transiting 
systems, providing additional radial velocity measurements 
for 10 systems described above, and include 10 systems with- 
out further reanalysis of orbital ephemeris. These systems 
are listed in Table |19| In Section |4.5[ we include a further 
16 systems, most of which had been discovered in the mean 
time. The planets involved in this study are listed on the 
webpage http://www.inscience.ch/transits/, where we also 
include the parameters usini (the projected rotation veloc- 
ity of the host star), Prot (the orbital period of the host star) 
and the projected spin-orbit angle A where available. 



Given the rapid rate of announcement of new transit- 
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Table 3. HARPS radial velocity measurements for WASP-2 (er- 
rors include random component only). 
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Table 4. HARPS radial velocity measurements for WASP-4 (er- 
rors include random component only). 
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4768, 
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Table 5. HARPS radial velocity measurements for WASP-5 (er- 
rors include random component only). 
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-29, 
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-29, 


,4994 
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-29, 
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-29, 


,3948 
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,5636 
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-29, 
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Table 6. HARPS radial velocity measurements for WASP-7 (er- 
rors include random component only). 



Time 
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[BJD-2450000] 


[km 


s-i] 


[km s-l] 


5003.41234 


-1 


.3253 


0.0206 


5005.47636 


-1 


..3669 


0.0122 


5006.50185 


-1 


.3244 


0.0123 


5007.42327 


-1 


.4822 


0.0129 


5008.39084 


-1 


.4143 


0.0131 


5009.38832 


-1, 


.3465 


0.0132 


5010.39331 


-1 


.4711 


0.0131 


5011.42884 


-1 


.4181 


0.0129 


5012.46735 


-1 


.3504 


0.0131 


5013.45923 


-1 


.4487 


0.0128 


5014.43881 


-1 


.4106 


0.0123 


5015.48148 


-1 


.3212 


0.0124 


5016.41444 


-1 


.4666 


0.0119 



Table 7. SOPHIE Radial velocity measurements for HAT-P-4 
(uncertainties include random component only). 
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[km 


s-1] 


[km 


s-1] 


5002.48517 


-10, 


,2995 


0, 


.0100 


5003.52118 


-10, 


,6910 


0, 


.0103 


5004.59910 


-10, 


,2681 


0, 


.0137 


5005.49926 


-10, 


,6377 


0, 


.0101 


5006.55335 


-10, 


,2564 


0, 
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-10, 
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5008.47624 


-10, 


,4027 





.0106 


5010.43095 


-10, 


,5681 


0, 


.0102 


5011.52259 


-10, 


,3835 


0, 


.0102 


5013.60648 


-10, 


.3090 


0, 


.0093 


5014.57426 


-10, 


,6862 


0, 


.0101 


5015.58518 


-10, 


,2680 


0, 


.0103 


5016.54123 


-10, 


,6808 





.0084 



Table 8. SOPHIE radial velocity measurements for HAT-P-7 
(uncertainties include random component only). 
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Time RV o-rv 

[BJD-2450000] [km s-i] [km s'l] 



5005.57207 


-0, 


.4489 


0.0107 


5006.57091 


-0, 


.2231 


0.0100 


5007.57482 


-0, 


.2619 


0.0105 


5008.44948 


-0, 


.4742 


0.0106 


5010.44234 


-0, 


.4505 


0.0107 


5011.51233 


-0, 


.2499 


0.0110 


5013.59448 


-0, 


.4207 


0.0114 


5014.56301 


-0, 


.1604 


0.0108 


5015.57356 


-0, 


.5090 


0.0107 


5016.55442 


-0, 


.2197 


0.0110 



Table 9. SOPHIE radial velocity measurements for TrES-2 (un- 
certainties include random component only). 



Time 


RV 




[BJD-2450000] 


[km s~^] 


[km s^"*"] 


4878.41245 


46.7905 


0.0091 


4879.38681 


46.9667 


0.0084 


4886.39349 


46.7748 


0.0095 


4887.44867 


46.9583 


0.0084 


4888.47514 


46.7722 


0.0085 


4889.40965 


46.8778 


0.0086 


4890.46546 


46.8994 


0.0085 


4893.41643 


46.8202 


0.0087 


4894.44335 


46.8073 


0.0121 



Table 10. SOPHIE radial velocity measurements for XO-2 (un- 
certainties include random component only). 
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3 ANALYSIS 

We used the radial velocity data , as well as the constraints 
on the orbital period P and mid-transit time Ttr (and e cos uu 
where available) from photometry as described in Section [5] 
To calculate the median values of the derived parameters 
and their corresponding uncertainties, we marginalise over 
their joint probability distribution using a Markov Chain 
Monte Carlo analysis with the Metropolis-Hastings algo- 
rithm. This has been described in the past by|Holman et al 



( 2006[) and our implementation is described in 



Pont et al. 



(2009b I. We model the radial velocity using a Keplerian or- 



bit and run the MCMC for 500,000 steps, the first 50,000 of 
which are then dropped to allow the MCMC to lose memory 
of the initial parameters. We verify that the autocorrelation 
length of each chain is much shorter than the chain length to 
that ensure the relevant region of parameter space is prop- 
erly explored. 

Although this procedure is common practice in the com- 
munity, we bring two changes. The first is a modification to 
the merit function that is used to work out the likelihood of 
a set of parameters given the data, to include the effects of 
correlated noise. This is described in Section [3.11 The sec- 
ond modification we bring is that we consider not only the 
case of an orbital model with a free eccentricity e, but we 
also work out the likelihood for a circular orbit (i.e. with e 
fixed at zero). We then compare the two models, by includ- 
ing a penalty for the additional complexity in the eccentric 
one (i.e. two additional degrees of freedom). We do this by 
using the Bayesian Information Criterion, as described in 
Section 13.31 

We report the median value in the chain for each pa- 
rameter, as well as the central 68.4% confidence interval on 
the parameter. A circular orbit model for an orbit that is 
in fact eccentric would artificially make the uncertainties in 
the derived parameters smaller, so in the case of the sys- 
temic velocity Vb and the semi-amplitude K, we report the 
median values from a circular orbit model, yet we include 
the confidence intervals derived from the eccentric model. 



3.1 The treatment of correlated noise 

Correlated noise can be important in the analysis of tran- 



sit light curves (Pont et al. 20061, and we included this in 



the analysis of radial velocity measurements in the case of 



WASP-12 (see Husnoo et al. 20111. If we assume uncorre- 
lated Gaussian noise when analysing data that is affected by 
correlated noise, we run the risk of overestimating the im- 
portance of a series of measurements that were obtained in 
quick succession, and this can have implications for example 
in estimating the orbital eccentricity. 

From Sivia (20061, the likelihood function for some 



data, given a model, is given by: 



P{-D\e,i) 



exp[-|(F-D)^C-^(F-D)] 
^(27r)^det(C) 



(3) 



where D is the radial velocity time series data expressed 
as a vector, 9 is the vector of model parameters, F is the 
predicted values from the Keplerian model. In this case, 
is defined by 



where C is the covariance matrix, which remains constant 
throughout the MCMC analysis for each system. In the case 
of independent measurements, the components of the covari- 
ance matrix C would be obtained using 



Cfe 







for k — k' 
otherwise. 



(5) 



whereas in the presence of some correlated noise, we modify 
this to include a squared exponential covariance kernel so 
that 



= 5k,k'(^k + ^ (^i exp ^ 



It} 



(6) 



where Gk is the formal uncertainty on each measurement k 
as obtained from the data reduction for that measurement, 



and the sum over M terms having the form g\ exp ^ 

allows us to include a number of stationary covariance func- 
tions to account for correlations in the noise, occuring over 
the timescales of hours to days. 

In practise, it can be tricky to estimate the values of 
Ti and Gi for radial velocity datasets, especially where the 
number of measurements is few or the phase-coverage is in- 
complete. Given the small datasets, we elect to use a sin- 
gle time-invariant correlation term, setting M = 1. The 
timescale is now called r, and the corresponding value of 
Oi is now called (Jr, where the subscript r indicates "red 
noise" (Pont et al. 20061. A fully Bayesian analysis would 



require that we assign priors to these two parameters and 
then marginalise over them. In practise, the sparse sampling 
and small datasets for radial velocity observations mean that 
it is very difficult to perform Bayesian marginalisation over 
these two parameters, and the results would depend on the 
prior space chosen (eg: r, logr, etc). We therefore use a sin- 
gle pair of parameters for r and Or- Using a single term in 
the sum in Equation [6] makes the expression less flexible, 
but prior experience shows that correlations over the ~ 1 
day timescale are particularly important for radial veloc- 
ity datasets, especially for measurements taken in the same 
night (see for example, Husnoo et al.|[20TT 



X 



D)^C"'(F- 



D) 



(4) 



for the case of 

WASP-12). For datasets where the reduced for a given 
model (circular or eccentric) was larger than unity, we es- 
timated Ot by repeating the MCMC analysis with different 
values of Or until the the best-fit orbit resulted in a reduced 
of unity for some optimal value of Gr- We discuss the 
estimation of r in the next sub-section. 



3.2 Estimation of r 

There is a degeneracy between Gr and r for the time sam- 
pling typical of our RV data: if we assume a long timescale 
compared to the interval of time between the measurements, 
we are asserting that we have a reason to believe that 
several measurements may have been systematically offset 
in the same direction. A measurement that occurs within 
that timescale but is offset to a very different extent from 
nearby measurements (e.g. if the correlation timescale r has 
been overestimated) will require a larger value of g^ for the 
dataset as a whole to yield a reduced of unity. 

To estimate r, we looked at several datasets for each 
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of the instruments HARPS, HIRES and SOPHIE. We re- 
peated the analysis in Section |3.1| using values of r in the 
range 0.1-5 d, to check for weather-related correlations. To 
see the effects of choosing between an eccentric orbit or a cir- 
cular orbit on our estimation of r, we carried each analysis 



twice, by adjusting cjr (see Section 3.1 1 to obtain a reduced 



X of unity (within 0.5%) for each orbital model (circular 
and eccentric). We plotted the optimal values of Or against 
T for several objects using data obtained from different in- 
struments separately, as shown in Figure [T] For WASP-2, 
we used our new HARPS me asurements (Table [3]) as well as 
SOPHIE measurements from [Collier Cameron et al. (20071. 
For WASP-4 and WASP-5, we used our new HARPS mea- 
surements (Tables |4] and [5|, and for HAT-P-7, we used our 
new SOPHIE measurements (Tab le [s]) a s well as HIRES 



measurements from [Winn et al. (2009c I. We found that 



for those datasets and objects where the orbital elements 
were well-constrained the plot showed a gentle increase in 
iTr with r, for r ^ 1.5d, then increased much faster for these 
datasets at a timescale of r > 1.5 d. For objects that have 
been observed with multiple instruments, this characteristic 
timescale is independent, both of the instrument used or the 
assumption about the eccentricity (i.e. free eccentricity or e 
fixed at zero), suggesting that the correlated noise is prob- 
ably related to weather conditions. We therefore assumed a 
correlation timescale of 1.5 d in the rest of this study, unless 
otherwise noted. This means that we are accounting for the 
red noise in the same-night measurements, and for measure- 
ments that are taken further apart in time, this procedure 
reduces to the more familiar "jitter" term. The value of ar 
inferred at r = 1.5 d in some cases varies by a few percent 
depending on the model chosen, i.e. eccentric or circular, 
and varies across datasets, as discussed later. 

We also investigated the effects of varying r on our final 
results. For the same systems discussed above, we plotted 
the 95% upper limit on the eccentricity as obtained from 
each dataset separately. The results are shown in Figure [2j 
where it is clear that the choice of r has no effect on the 
final result for r ^ 1.5 d. The only exception is WASP- 
2 (HARPS), where we only have 8 measurements and the 
phase coverage is not as complete as for the other objects 
(see Figure [To| ) . Similarly, the derived parameters Vo and K 
did not vary appreciably with r. 

3.3 Model selection 

Determining whether an orbit is consistent with circular 
(Model 1) or eccentric (Model 2), is an exercise in model 
selection. If we assume the prior probability of the circular 
and eccentric models are the same, we can use the Bayesian 



Information Criterion (BIC) ( Liddle 2007 1 to decide between 



the two models. This is equivalent to working out the Bayes 
factor P(dataiModeli)/P(data|Model2), subject to the as- 
sumptions described below. The Bayes factor is the ratio of 
marginal likelihoods for each model, each of which is given 
from 

P(data|Modelj) = / L{Qj\A&ta) y. P{ej\MoAe\j)dQj (7) 

where Oj represent the vector of parameters for each model 
j, Z/(0-,|data) is the likelihood and P(Oj|Modelj) is the joint 
posterior distribution of the parameters. 



As described in Section [3] above and in I Pont et al.l 
( |2009b| l, the MCMC process produces the joint posterior 
distribution for the parameters, and we also obtain a max- 
imum likelihood Lmax, corresponding to the smallest value 
of (as given in Section 3.1 1 for each model. We then use 



the Bayesian Information Criterion (Liddle 20071 as given 
by. 



BIC = -21nLn,ax + fclniV, 



(8) 



where N is the number of measurements, k is the number of 
parameters in the model used. This simplifies the expression 
for the marginal likelihood by performing the integration us- 
ing Laplace's method and assumes a flat prior. If we replace 
imax with the expression given by P(D|S,/) in equation [s] 
above, 



BIC = Xmin + fclnA^ + ln((27r)^|C|' 



(9) 



where Xmin is the minimum value of x^ achieved by the 
model, A'^ is the number of measurements, k is the number 
of parameters in the model, and |C| is the determinant of 
the correlation matrix given in Section [3. 1| above. 

The radial velocity data for a Keplerian orbit involves 6 
free parameters: the period P, a reference time such as the 
mid-transit time Ttr, a semi-amplitude K, a mean velocity 
offset Vb, the argument of periastron u) and the eccentricity 
e. Following [Fbrg| ( [20061 ), we use the two projected compo- 
ments ecosw and esinoj instead of e and uj, to improve the 
efficiency of the MCMC exploration. In this study, we use 
the period P, mid-transit time Ttr and their corresponding 
uncertainties as a priori information. We thus count them as 
two additional measurements in the calculation of the BIC, 
while the number of free parameters in each model (circular 
or eccentric) is now decreased by two. In this circular 
model would have 2 free parameters (Vb and K), while an 
eccentric model would have 4 free parameters (Vb, K, e, and 
a;). 

The term k\nN thus penalises a model with a larger 
number of parameters (for example, an eccentric orbit), and 
we seek the model with smallest BIC. For each object, we 
repeated the MCMC analysis using the optimal value for 
cjr for a circular orbit and an eccentric orbit separately, at 
r = 1.5 d, unless otherwise noted. We call these two families. 
For each family, we performed a fit with a circular model and 
an eccentric model. In most cases, the two families agreed 
on a circular model (indicated by "C" in Table |19[ ) or an 
eccentric model (indicated by "E" if e > 0.1 in Table 191, 



indicating this with a smaller BICc or a smaller BICe respec- 
tively. If the two families favoured a circular (or eccentric) 
orbit, we give the parameters from the family using an op- 
timal value of Gr for the circular (or eccentric) orbit. In a 
number of such cases, however, the upper limits on the or- 
bital eccentricity were larger than e = 0.1. We labelled these 
eccentricities as "poorly constrained" (indicated by "P" in 
Table [Tol l. In a few cases, the small number of measurements 
or the quality of measurements (e.g. for faint targets, or low 
mass planets) meant the two families disagreed: the fam- 
ily using the optimal value of (Jr for a circular orbit gave a 
smaller value of BICc, favouring the circular orbit and the 
family using the optimal value of ar for an eccentric orbit 
gave a smaller value of BICe, favouring the eccentric orbit. 
We labelled these cases "poorly constrained" as well. 
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Figure 1. Plot showing the degeneracy between ar and t for objects where the orbital parameters are well constrained (see Sections |3.1| 
and |3.2[ l because a sufficient number of radial velocity measurements is available and provides sufficient phase coverage. The object being 
studied is shown in the title for each panel, and the instrument used for the measurements are shown in parenthesis. As can be seen on 
each plot, the optimal ar that gives a reduced of unity for each dataset increases slowly with t for r < 1.5 d, but increases faster 
after 1.5 d. This hints that the systematic effects occur on a timescale of 1.5 d, and could be related to the weather. Note the SOPHIE 
data for WASP-2 did not provide full phase coverage — the solution did not converge for an eccentric model and the knee at r ~ 1.5 d 
is less pronounced. 



4 RESULTS 

The results of this study are shown in Table |19[ We place 
constraints on the eccentricities of transiting planets for 
which enough data is available. We analysed radial veloc- 
ity data for 54 systems. For 8 systems, we used our new 
radial velocity data (described in Section [2| in addition to 
existing RVs from the literature. For the other 46 systems, 
we reanalysed existing RVs from the literature. 

In Section [4. 1[ we describe the planets for which we do 
not consider the evidence for an orbital eccentricity com- 
pelling, despite previous evidence of a departure from circu- 



larity (e > 1(7 from zero), followed by Section [4.2[ where we 
describe the planets for which we consider the orbital eccen- 
tricity to be either so small as to be undetectable or com- 
patible with zero. In Section |4.3| we describe planets that 
can be safely considered to be on eccentric orbits and fi- 
nally, in Section |4.4| we describe the planets for which we 
consider the orbital eccentricity to be poorly constrained (as 
described in Section 3.31. In the following Sections we also 



include a discussion of the evidence for eccentricity for 26 
other systems from the literature. 
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Figure 2. Plot showing tlie effect of varying r on the 95% upper limit on the derived eccentricity for WASP-2, HAT-P-7, WASP-4 and 
WASP-5. Except for WASP-2, where the phase coverage of the HARPS data is incomplete (see Figure [To] l, varying r has no effect on the 
95% upper limit of the derived eccentricity for timescales of a few days, r > 1.5 d. Note: each line on this plot is made a single dataset. 



4.1 Planets with orbits that no longer qualify as 
eccentric according to this study 

In a number of cases in the past, the derived eccentricity 
from an MCMC analysis deviated from zero by more than 
Icr, for example CoRoT-5b, GJ436b, WASP-5b, WASP- 
6b, WASP-lOb, WASP-12b, WASP-14b, WASP-17b and 
WASP-18b. In this Section, we discuss the cases of 7 planets, 
CoRoT-5b, WASP-6b, WASP-lOb, WASP-12b, WASP-17b, 
WASP-18b and WASP-5b, that are shown to have orbital 
eccentricities that are compatible with zero. 

CoRoT-5 

CoRoT-5b is a 0.46 Mj planet on a 4.03 day orbit around 
a F9 star (V=14.0), first reported by [Rauer et al.| (|2009 |. 



Using 6 SOPHIE measurements (one of which is during the 
spectroscopic transit, which we ignore in this study) and 
13 HARPS measurements, the authors derived a value of 
eccentricity e = 0.09to°ol- 

In our study, we used the formal 
uncertainties quoted with the data without any additional 
noise treatment, since they resulted in a reduced less than 
unity for both an eccentric and a circular orbit. We imposed 
the prior information from photometry P — 4.0378962(19) 
and Ttr = 2454400.19885(2) from the [Rauer et al.| ( |2009 l 
and obtained a value of Xc = 15.97 for the circular orbit and 
a value of Xe = 13.50 for the eccentric orbit (e = 0.086tQ'.o545 
e < 0.26). Using TV — 20, k — 3 and fc = 5 for the circular 



(two datasets, each with one Vb and a single K) and eccentric 
orbits respectively, we obtained BICc = 151.05 and BICe ~ 
154.57. A smaller value of BICc means the circular orbit 
cannot be excluded. 

WASP-6 

WASP-6b is a 0.50 Mj planet on a 3.36 day orbit around 
a G8 star (V=11.9), first reported by |Gillon eTel] ^009\ . 
Using 35 CORALIE measurements and 44 HARPS mea- 
surements (38 of which occur near or during a spectroscopic 
transit, which we ignore in this study), the authors derived a 
value of eccentricity e — 0.054:t.Q'Q\^- In our study, we used 
the 35 CORALIE measurements and the 6 HARPS mea- 
surements that were not taken in the single night where the 
spectroscopic transit was observed. We used (Jr = m s~^ 
for CORALIE (the data produces a reduced x^ = 0.89 when 
fitted with a circular orbit, indicating overfitting) but for 
HARPS we used r = 1.5 d and ar = 4.15 m s~^ to obtain 
a reduced x^ of unity for the circular orbit. We obtained 
a value of Xc = 38.09 for the circular orbit and a value 
of Xe = 33.58 for the eccentric orbit (e = 0.041 ± 0.019, 
e < 0.075). Using A'" = 43 (41 RVs and two constraints from 
photometry), fe = 3 and fe = 5 for the circular (two datasets, 
each with one Vo and a single K) and eccentric orbits re- 
spectively, we obtained BICc = 333.25 and BICc ~ 336.27. 
We repeated the calculations, using ar = for CORALIE 
(the data produces a reduced x^ = 0.85 when fitted with an 
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eccentric orbit, indicating overfitting) but for HARPS we 
used r = 1.5 d and ct^ = 3.59 m to obtain a reduced 
of unity for the eccentric orbit. We obtained a value of 
Xc = 39.20 for the circular orbit and a value of Xe ~ 34.47 
for the eccentric orbit (e = 0.043 ± 0.019, e < 0.075). Using 
N = 43, k = 3 and fc = 5 for the circular and eccentric orbits 
respectively, we obtained BICc ~ 333.60 and BlCe = 336.39. 
We therefore find that the circular orbital solution cannot 
be excluded, but the possibility that e > 0.1 is rejected. 

WASP-10 

WASP-lOb is a 2.96 Mj planet on a 3.09 day orbit around a 



clear to us whether or not the orbital eccentricity is non-zero 
as claimed in 



Christian et al 



K5 star (V=12.7), first reported by Christian et al. ( 2009 1 



Using 7 SOPHIE measurements and 7 FIES measurements, 
the authors derived a value of eccentricity e = 0.059to oo4- 
The FIES data yielded a reduced x^ l^ss than unity with 
both eccentric and circular orbits, indicating overfitting, so 
we set (Tr = m s~^. 

For the SOPHIE data, used r = 1.5 d, cr^ = 54.5 m s"^ 
to obtain a reduced x^ of unity for the circular orbit. We 
reanalysed all the radial velocity measurements, and applied 
the prior from photometry P = 3.0927636(200) and Ttr = 
2454357.8581(4) from .Christian et al.^ (2009) . We obtained 
a value of Xc = 13.49 for the circular orbit and a value of 
Xe = 7.47 for the eccentric orbit (e = 0.049 ± 0.022, less 
significant than the original claim). Using 14 measurements 
and two priors from photometry (A'^ — 16), k — 3 and fc = 5 
for the circular (two datasets, each with one Vb and a single 
K) and eccentric orbits respectively, we obtained BICc ~ 
151.50 and BlCe = 151.01. This now appears to show only 
a marginal support for an eccentric orbit. 

We plotted the SOPHIE radial velocity data against 
time, as shown in Figure [3] and overplotted a circular or- 
bit as well as an eccentric orbit. Due to the long time be- 
tween the first two measurements and the last five, we plot 
them in separate panels, shown on the left and right respec- 
tively. It is clear that the first measurement is pulling the 
eccentricity upwards, and we suspect from experience that 
the long term drifts in the SOPHIE zero point in HE mode 
for faint targets could have affected the first two measure- 
ments. We therefore repeated our calculations using only 
the last five measurements from the SOPHIE dataset and 
the whole FIES dataset, and set ar = 45.5 m s~^ for SO- 
PHIE. This time, we obtained a value of Xc = 11-81 for 
the circular orbit and a value of Xe ~ 7.64 for the eccentric 
orbit (e — 0.043 ± 0.035). Using 12 measurements and two 
priors from photometry (TV = 14), fc = 3 and fc = 5 for 
the circular (two datasets, each with one Vb and a single K) 
and eccentric orbits respectively, we obtained BICc ~ 128.71 
and BICc = 129.83, this time favouring the circular orbit. 
We repeated this calculation, and set (Tr = m s^^ for 
both SOPHIE and FIES, as each dataset gave a reduced x^ 
of less than unity for the eccentric orbit (SOPHIE reduced 
X^ = 0.64, FIES reduced x^ = 0.45). This time, we obtained 
a value of Xc = 19.30 for the circular orbit and a value of 
Xe = 10.65 for the eccentric orbit (e = 0.080 ± 0.055). Using 
12 measurements and two priors from photometry (A'^ = 14), 
fc = 3 and fc = 5 for the circular (two datasets, each with 
one Vb and a single K) and eccentric orbits respectively, 
we obtained BICc = 128.33 and BICe = 124.97, this time 
favouring the eccentric orbit once again. It is therefore un- 



(20091 



'Maciejewski et al. (20111, used transit timing varia- 
tion analysis and reanalysed the radial velocity data, to 
obtain an eccentricity that is indistinguishable from zero 
(e — 0.013 ± 0.063). They argued instead that the original 
detection of an eccentricity had been influenced by starspots. 
The difference between our value of eccentricity and that 



derived by Maciejewski et al. (20111 is probably due to the 



fact that the latter used a two planet model, which can re- 
duce the derived eccentricity further — sparse sampling of 
the radial velocity from a two planet system can lead to an 
overestimated eccentricity. 

WASP- 12 

WASP- 12b is a 1.41 Mj planet on a 1.09 day orbit around 



a F9 star (V=11.7), first reported by Hebb et al. (20091 



Using SOPHIE measurements, the original authors derived 
a value of eccentricity e = 0.049±0.015. [Husnoo et al.| ( |20lT| 
used new SOPHIE radial velocity measurements, as well as 



the original transit photometry from Hebb et al. ( 2009 1 and 
the secondary eclipse photometry from |Campo et al.| ( |201l] ) 
to suggest that the eccentricity was in fact compatible with 
zero (e = 0.017lHlo)- 

WASP- 17 

WASP-17b is a 0.50 Mj planet on a 3.74 day orbit around a 
F6 star (V=11.6), first reported by [Anderson et al.| ( [20Tol ). 
Using 41 CORALIE measurements (three of which are dur- 
ing the spectroscopic transit, which we ignore in this study) 
and 3 HARPS measurements, the authors considered three 
cases: first imposing a prior on the mass M* of the host 
star, secondly imposing a main-sequence prior on the stellar 
parameters and thirdly with a circular orbit. They derived 
values of eccentricity e = 0.129115:^68 ^nd e = 0.237j:!5;!J^^ 
for the first two cases respectively. We set ar = for both 
HARPS and CORALIE since we obtained a reduced x^ of 
slightly less than unity for both eccentric and circular orbits 
for either dataset alone, indicating overfitting. We obtained 
a value of Xc = 37.98 for the circular orbit and a value 
of Xe ~ 35.94 for the eccentric orbit. Using 41 measure- 
ments and two priors from photometry {N — 43), fc = 3 and 
fc = 5 for the circular (two datasets, each with one Vb and 
a single K) and eccentric orbits respectively, we obtained 
BICc = 399.31 and BICc = 404.80. We thus find that the 
circular orbit cannot be excluded, agreeing with the third 



case (e — 0, fixed) considered in Anderson et al. (20101 and 
rejecting the two derived values of eccentricity in that paper. 

WASP-18 

WASP-18b is a 10.3 Mj planet on a 0.94 day orbit ar ound a 
F6 star (V=9.3), first reported by [Hellier et~aL] ( |2009t . Using 
9 CORALIE measurements (we drop the third measurement 
in our final analysis, since it produces a 5-cr residual that is 
not improved by an eccentric orbit, suggesting that it is a 
genuine outlier), the authors derived a value of eccentricity 
e = 0.0092 ± 0.0028. In our study, we set r = 1.5 d and 
ar — 20.15 m s~^ to obtain a reduced x^ of unity for the 
circular orbit. We obtained a value of Xc ~ 8.17 for the 
circular orbit and a value of Xe ~ 6.64 for the eccentric orbit 
(e = 0.007 ± 0.005, e < 0.018). Using TV = 10, fc = 2 and 
fc = 4 for the circular (one dataset, with one Vb and a single 
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Figure 3. Plot showing SOPHIE radial velocity data from |Christian et al.| | |200"9l l for WASP-10, plotted against time. The plot has been 
split along the time axis into two panels (left and right) to remove the 160 days without measurements, for clarity. A circular orbit 
(solid line) and an orbit with the best-fit eccentricity (e = 0.048) are overplotted. The residuals relative to the circular orbit are shown 
in the bottom panel. 



K) and eccentric orbits respectively, we obtained BICc ~ 
75.34 and BICe = 78.41. We repeated the calculations using 
CTj. = 22.5 m s~^ to obtain a reduced of unity for the 
eccentric orbit. We obtained a value of Xc ~ 7.14 for the 
circular orbit and a value of Xe ~ 6.00 for the eccentric 
orbit (e = 0.008 ± 0.005, e < 0.019). Using iV = 10, A: = 2 
and fc = 4 for the circular (one dataset, with one Vb and 
a single K) and eccentric orbits respectively, we obtained 
BICc = 75.50 and BIC^ = 78.97. We thus find that the 
circular orbit cannot be excluded, in contrast to |Hellier et al.| 
(20091). The possibility that e > 0.1 is excluded. 



WASP-5 (new HARPS data) 

WASP-5b is a 1.6 planet on a 1.63 day orbit around a 



G4 star (V=12.3), first reported by Anderson et al. ( 2008 1 



Gillon et al. (2009b I used z-band transit photometry from 



the VLT to refine the eccentricity to e = 0.038lQ g^|, and 
the authors made a tentative claim for the detection of a 
small eccentricity. We analysed our 11 new HARPS mea- 
surements for WASP-5 and the 11 CORALIE RVs from UiTl 



derson et al. 



( 2008 \ using the photometric constraints on 
1.6284246(13) and mid-transit time 



the orbital period P 
Ttr = 2454375.624956(24) from |Southworth et al.| ( [20()9 l 
We use r 



= 1.5 d. 



10.6 m s"' for HARPS and 



(Tr = 4.3 m s~ for CORALIE to obtain a value of reduced 
X^ of unity for the circular orbit for each dataset separately. 
We ran the MCMC twice: the first time fitting for the sys- 
temic velocity no and semi-amplitude K, and the second 
time adding two parameters ecosoi and esint^; to allow for 
an eccentric orbit. The best fit result is shown in Figure [4] 
The residuals for a circular orbit are plotted, and a signal is 
clearly present in the residuals. The value of x^ for the cir- 
cular orbit is 24.36 and that for an eccentric orbit is 20.57. 
This results in a value of BICc = 169.40 for the circular 
orbit and BICe ~ 171.97 for the eccentric orbit, given 22 
measurements, 2 constraints from photometry and 3 and 5 
free parameters respectively for each model. 

We repeated the above analysis using r = 1.5 d, 



(Jr = 9.4 m s~ for the HARPS dataset to obtain a value 
of reduced x^ of unity for the eccentric orbit and ar — 
m s~^ for CORALIE (which resulted in a reduced x^ of 
0.58). This time, we obtained a value of x^ for the circu- 
lar orbit is 27.35 and that for an eccentric orbit is 23.00. 
This leads to a value of BICc ~ 170.37 for the circular orbit 
and BICe = 172.38 for the eccentric orbit. Once again, the 
circular orbit is favoured. 

A keplerian model, circular or eccentric (e = 0.012 ± 
0.007) does not account for the scatter in the data the 
HARPS dataset as shown in Figure [4] We have therefore 
plotted the radial velocity measurements, the bisector span, 
the signal to noise at order 49, the contrast and full width 
at half maximum for the cross-correlation function against 
the same time axis. The trend in radial velocity residuals 
can be seen to be correlated with both the bisector span 
and the full width at half maximum of the cross correlation 
function. This suggests a line shape change that's related 
to either weather effects or instrumental systematics. The 
timescale of this variation is compatible with both scenar- 
ios. The bisector inverse span is generally directly correlated 
with the residuals, which weighs against a scenario involv- 
ing stellar activity, but this is not so clear for the first three 
measurements — the drift could be due to stellar activity or 
an additional planetary or stellar companion. 

We extended the model with a linear acceleration of the 
form 



V{t) = t^kcplorian (t) +j{t-to), 



(10) 



and fitted the HARPS data alone using to = 2454768 
(to allow the MCMC to explore values of -y more efficiently) 
and reran the MCMC twice: once for a circular orbit and 
once for an eccentric orbit. Firstly, we used ar = 10.6 m s~^ 
for the HARPS dataset, and the linear trend for a circular 
orbit resulted in 7 = —2.6 ± 2.9 m s~^ yr^^ and that for an 
eccentric orbit is 7 = —2.0 ± 2.9 m s~^ yr~^. The best fit 
result is shown in Figure |5] and the residuals for a circular 
orbit are plotted in the bottom panel. The value of for the 
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Figure 4. HARPS measurements of WASP-5 plotted against time (left) and phase with respect to the mid-transit time Ttr (right). 
In each case, a soUd Une is overplotted to represent a circular orbit and the residuals are plotted for this circular orbit. It is clear that 
a signal is present in the residuals (see text). An eccentric orbit with the best-fit value of e = 0.012 is overplotted in both panels with a 
dotted line, but it is indistinguishable from the circular solution at this scale. Note the trend that is apparent in the residuals (second 
panel from the top on both the time and phase plots). We correct for this using a linear acceleration term in our model (see Figure [sjl. 



circular orbit is 10.24 and that for an eccentric orbit is 7.70. 
This results in a value of BlCc.iin = 71.91 for the circular 
orbit and BICe, zin = 74.49 for the eccentric orbit, given 11 
(N=13) measurements, 2 constraints from photometry and 
3 and 5 free parameters respectively for each model. We 
repeated these calculations using ar — 9.4 m s~^ for the 
HARPS dataset, and the linear trend for a circular orbit 
resulted in 7 = —3.7 ± 1.3 m s~^ yr~^ and that for an 
eccentric orbit is 7 = —3.3 ± 1.3 m s~^ yr~^- The value of 

for the circular orbit is 15.46 and that for an eccentric 
orbit is 13.50. This leads to a value of BICc.Kn = 69.72 
for the circular orbit and 'SlCe,Hn = 72.89 for the eccentric 
orbit. The circular orbit is not excluded, and the possibility 
that e > 0.1 is excluded. The results for both models, one 
including the linear trend but excluding the CORALIE data, 
and one including the CORALIE data but excluding the 
linear trend are shown in Table |11| In both cases, we give 
results for the case where Gr is chosen to yield a reduced 

of unity for the circular orbit. We attempted to repeat 
this using both the CORALIE and HARPS datasets, but we 



were unable to obtain a fit with the MCMC, because of the 
long time scale between the two datasets. 



4.2 Planets on circular orbits 

We establish that 20 planets have orbital eccentricities com- 
patible with zero and the 95% upper limits are smaller than 
695 = 0.1. In this Section, we describe the planets WASP-4b, 
HAT-P-7b, TrES-2 and WASP-2b, for which we introduce 
new RVs. We also establish that the 95% upper limits on 
the eccentricities of WASP-5b, WASP-12b and WASP-18b, 
which have been described in Section 14.11 above. In addi- 
tion, we give the 95% upper limits on the eccentricities of 
CoRoT-lb, CoRoT-3b, HAT-P-8b, WASP-3b, WASP-16b, 
WASP-19b, WASP-22b, WASP-26b and X0-5b in Table[l9] 
We discuss the evidence for circular orbits for HAT-P-13b, 
HD 189733b, HD209458b and Kepler-5b at the end of this 
section. 

WASP-4 (new HARPS data) 

WASP-4b is a 1.2 Mj planet on a 1.34 day orbit around a 
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Figure 5. HARPS measurements of WASP-5 plotted against time (left) and phase with respect to the mid-transit time Ttr (right). 
In each case, a solid line is overplotted to represent a circular orbit and the residuals are plotted for this circular orbit. It is clear that 
a signal is present in the residuals (see text). A model for an eccentric orbit with the best-fit value of e = 0.013 is overplotted in both 
panels with a dotted line, but it is indistinguishable from the circular solution at this scale. Both include the linear trend (see Section [4.2| 
and Figure [4]|. 



Parameter 



Anderson et al. 



I 2008 1 



HARPS only, this work 
(with linear trend) 



HARPS & CORALIE, this work 
(no linear trend) 



Centre-of-mass velocity Vq [m s~^] 
Orbital eccentricity e 
Argument of periastron u [°] 
e cos ui 
e sintj 

Velocity semi-amplitude K [ms~^] 



20010. 5±3.4 
(adopted) 
(unconstrained) 



277.8±7.8 



20018±12 
0.013 ±0.008 (< 0.029) 
(unconstrained) 
0.002±0.003 
0.012±0.010 
266.4±1.3 



20009. 9±7.4 (HARPS) 
0.012 ±0.007 (< 0.026) 
(unconstrained) 
0.003±0.003 
0.011±0.009 
266.9±1.3 



Table 11. System parameters for WASP-5. Left: [Anderson et al.| | [2008[ l. Right: Results from our HARPS radial velocity data alone, and 
results from using both our HARPS data and the original CORALIE data in [Anderson et al.| ( |2008| l. Median values for Vb and K are 
quoted for the circular orbits, as well as 68.3% confidence limits obtained from the eccentric solution (see section Analysis). 



G7 star (V=12.5), first reported by Wilson et al. ( 2008 1 



We analysed our 14 new HARPS measurements and the 



14 CORALIE measurements from Wilson et al. ( 2008 1 for 



WASP-4 and used the photometric constraints on the or- 
bital period P — 1.33823214(71) and mid-transit time Ttr = 
2454697.797562(43) from [Winn eTaT] ( [SOOg] ) . 

We estimate r = 1.5 d, cr^ = 11 m s"^ for the HARPS 
dataset and ar — 4.5 m s"'^ for the CORALIE dataset to ob- 
tain a reduced of unity for a circular orbit for each dataset 
separately. We ran the MCMC twice: the first time fitting 
for the systemic velocity vo and semi-amplitude K only, ie. 
a circular orbit {k = 2), and the second time adding two 
parameters ecoscj and esinoj to allow for an eccentric orbit 
(fe = 4). The best fit result is shown in Figure [6] The resid- 
uals for a circular orbit are plotted, and a signal is clearly 
present in the residuals. The value of for the circular or- 
bit is 27.13 and that for an eccentric orbit is 24.32. This 
leads to a value of BICc = 208.62 for the circular orbit 
and BICe ~ 212.55 for the eccentric orbit, given 14 mea- 



surements, 2 constraints from photometry and 2 and 4 free 
parameters respectively for each model. 

We repeated the calculations, estimating r — 1.5 d, 
(Jr = 10.1 m s~^ for the HARPS dataset and ar — 7.1 m s~^ 
for the CORALIE dataset to obtain a reduced of unity 
for an eccentric orbit for each dataset separately. The value 
of x^ for the circular orbit is 27.29 and that for an eccentric 
orbit is 24.33. This leads to a value of BICc = 208.72 for the 
circular orbit and BIC^ — 212.51. 

Note the trend that is apparent in the residuals in Fig- 
ure [6] We have therefore plotted the radial velocity mea- 
surements, the bisector span, the signal to noise at order 49, 
the contrast and full width at half maximum for the cross- 
correlation function against the same time axis. For most 
measurements, the trend in radial velocity residuals can be 
seen to be correlated with both the bisector span and the 
full width at half maximum of the cross correlation function. 
This suggests a line shape change that's related to either 
stellar activity, weather effects or instrumental systematics. 
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Figure 6. HARPS measurements of WASP-4 plotted against time (left) and phase with respect to the mid-transit time Ttr (right). 
In each case, a solid line is overplotted to represent a circular orbit and the residuals are plotted for this circular orbit. It is clear that 
a signal is present in the residuals (see text). An eccentric orbit with the best-fit value of e = 0.005 is overplotted in both panels with a 
dotted line, but it is indistinguishable from the circular solution at this scale. Note the trend that is apparent in the residuals (second 
panel from the top on both images). We attempt to correct for this by repeating our calculations with a linear acceleration term in the 
model (see Section |4.2| and Figure [vjl. 



The timescale of this variation is compatible with all three 
scenarios. 

We repeated the calculations for the HARPS dataset 
alone, and added a linear component to the radial velocity 
model in the same way we did for WASP-5 in Section |4.1| 
and we set to = 2454762 (to allow the MCMC to explore 
values of 7 more efficiently) and reran the MCMC twice: 
once for a circular orbit and once for an eccentric orbit. We 
set r = 1.5 d and err = 11 m s^^ for the HARPS dataset. 

The best fit result is shown in Figure [7] The residuals 
for a circular orbit are plotted, and a signal is clearly present 
in the residuals. The linear trend for a circular orbit results 
in 7 = 1023 ±490 m s"'^ yr~^ and that for an eccentric orbit 
is 7 = 919 ± 500 m s~^ yr"^ 

The value of for the circular orbit is 9.83 and that 
for an eccentric orbit is 7.51. This leads to a value of 



BICc = 92.02 for the circular orbit and BICe = 95.26 for 
the eccentric orbit, given 14 measurements, 2 constraints 
from photometry and 3 and 5 free parameters respectively 
for each model. 

We repeated the calculations, setting r = 1.5 d, ar = 
10.05 m s"^ for the HARPS dataset. The value of for 
the circular orbit is 10.30 and that for an eccentric orbit is 
8.07. This leads to a value of BICc = 91.18 for the circular 
orbit and BICe = 94.49. In all cases, the circular orbit is 
not excluded. The results for both models, one including 
the linear trend but excluding the CORALIE data, and one 
including the CORALIE data but excluding the linear trend 
are shown in Table [12] In both cases, we give results for the 
case where a,, is chosen to yield a reduced of unity for 
the circular orbit. We reject the possibility that e > 0.1. 
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Figure 7. HARPS measurements of WASP-4 plotted against time (left) and phase with respect to the mid-transit time Ttr (right). 
In each case, a soUd Une is overplotted to represent a circular orbit and the residuals are plotted for this circular orbit. It is clear that 
a signal is present in the residuals (see text). An eccentric orbit with the best-fit value of e = 0.004 is overplotted in both panels with a 
dotted line, but it is indistinguishable from the circular solution at this scale. The linear trend from Figure [6] has now been included in 
the model. 



Parameter 



Wilson et al. 1 2008 I 



HARPS only, this work 
(with linear trend) 



HARPS & CORALIE, this work 
(no linear trend) 



Centre-of-mass velocity Vq [m s^^] 
Orbital eccentricity e 
Argument of periastron w [°] 
ecosuj 
esino; 

Velocity semi-amplitude K [ms~^] 



57733±2 
(adopted) 
(unconstrained) 



240±10 



57773±10 
0.004 ±0.003 (<0.011) 
(unconstrained) 
0.004±0.003 
-0.002±0.004 
233.1±2.1 



57790.8±5.7 
0.005 ±0.003 (<0.011) 
(unconstrained) 
0.003±0.003 
-0.004±0.004 
233.7±2.0 



Table 12. System parameters for WASP-4. Left: [Wilson et al.| l |2008l l. Right: Results from our HARPS radial velocity data alone, and 
results from using both our HARPS data and the original CORALIE data in [Wilson et al."] ( |2008[ l. Median values for Vq and K axe 
quoted for the circular orbits, as well as 68.3% confidence limits obtained from the eccentric solution (see section Analysis). 



HAT-P-7 (new SOPHIE data) 

HAT-P-7b is a 1.8 Mj planet on a 2.20 day orbit around an 



F6 star (V=10.5), first reported by [Pal et al]p008| ). We use 
13 new SOPHIE radial velocity measurements and 16 out of 
the 17 HIRES measurements in Winn et al. t2009c) (we drop 
one in-transit measurement) to work out the orbital param- 
eters of HAT-P-7b. We impose th e period P = 2 . 20473 3(10) 
d as given from photometry in Welsh et al. (20101 and 
mid-transit time Ttr = 2454731 . 67929(4 3) BJD as given 

We set r 



= 1.5 
12.9 m s"^ for 



from photometry in Winn et al. (2009c 
d, ar = 9.41 m s"^ for HIRES and 
SOPHIE to obtain a reduced of unity for the best-fit cir- 
cular orbit for each dataset separately. We used 29 measure- 
ments in all, and count the two constraints from photometry 
as two additional data points to obtain A'' = 31, and used 
k — 4 for the circular orbit (two Vo, one for each dataset, 
the semi-amplitude K and a constant drift term 7, since 



Winn et al. ( 2009c I found evidence for a distant companion 
in the system and we set to = 2454342). We repeated this 
analysis with an eccentric orbit fe = 6 (4 degrees of free- 



dom for the circular orbit with two datasets and a Hnear 
acceleration, and 2 additional degrees of freedom for the ec- 
centricity, ecoso; and esintj). The orbital parameters are 
given in Table [Ts] and the radial velocity dataset is plotted 
in Figure [8] with residuals shown for a circular orbit. The 
Figure also shows models of a circular and an eccentric or- 
bit (with e = 0.014), but they are almost undistinguishable. 
For the circular orbit, we obtained = 26.94, and a value 
of BICc = 222.81 and for the eccentric orbit, we obtained 
= 23.98 and a value of BIC^ = 226.72. We repeated the 
calculations and set t = 1.5 d, (t^ = 8.2 m s~^ for HIRES 
and CT,. = 8.2 m s~^ for SOPHIE to obtain a reduced oi 
unity for the best-fit eccentric orbit. For the circular orbit, 
we obtained x^ = 35.65 and a value of BICc = 224.14 and 
for the eccentric orbit, we obtained x^ = 31.89 and a value 
of BICe = 227.25. We therefore find that the circular orbit 
cannot be excluded for HAT-P-7b. Further, we exclude the 
possibility that e > 0.1. 



TrES-2 (new SOPHIE data) 
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Parameter 



HIRES, Winn et al. (2009c 



HIRES+SOPHIE, this work 



Centre-of-mass velocity Vq [m s^-"^] 
Orbital eccentricity e 
Argument of periastron uj [°] 
e cos id 
esint*; 

Velocity semi-amplitude K [ms~^] 
Constant radial acceleration 7 [m 



-51.2±3.6 -49.96±6.0 (HIRES) and -10510±10 (SOPHIE) 

699% <0.039 0.014±0.010 (e < 0.038) 

— (unconstrained) 

-0.0019±0.0077 -0.007±0.004 

0.0037±0.0124 -0.011it0.015 

211.8±2.6 213.8±1.2 

21.5±2.6 21.1±4.2 



Table 13. System parameters for HAT-P-7. Left: | Winn et al.|p009c| . Right: Results from our SOPHIE radial velocity data. Median 
values for Vq a^iid K are quoted for the circular orbits, as well as 68.3% confidence limits obtained from the eccentric solution. The upper 
95% limit is also given for the eccentricity from our analysis. 
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Figure 8. Plot showing our new SOPHIE radial velocity data for HAT-P-7, plotted against time (left), and orbital phase (right) with 
respect to Ttr- A circular orbit (solid line) and an orbit with the best-fit eccentricity (dotted line, but almost undistinguishable from the 
circular solution since e = 0.014) are overplotted. The residuals relative to the circular orbit arc shown in the bottom panels. 



TrES-2b is a 1.3 Mj planet on a 2.47 day orbit around a GO 



star (V=11.4), first reported by O'Donovan et al. ( 2006 1 



We use 10 new SOPHIE radial velocity measurements and 



tiie 11 HIRES measurements in O'Donovan et al. (20061 



to work out the orbital parameters of TrES-2b. We impose 
the period P — 2.470614(1) d and mid-transit time Ttr = 



2453957.634 92(13) BJD as given from photometry in |Raetz 
et al.| ( |2009] ). 

We set r = 1.5 d and cr^ = 6.8 m s"^ for SOPHIE to 
obtain a reduced of unity for the best-fit circular orbit 
(using the SOPHIE data alone), and set ar — Om s~^ for 
the HIRES data since a circular orbit for that dataset alone 
yields a reduced of 0.72, indicating over-fitting. We used 
21 measurements in all, and count the two constraints from 
photometry as two additional datapoints to obtain A'^ = 23, 
and used A; = 3 for the circular orbit (two Vb, one for each 
dataset, and the semi-amplitude K). We repeated this anal- 
ysis with an eccentric orbit k = 5 (three degrees of freedom 
for the circular orbit, and two additional degrees of freedom 
for the eccentricity, ecosoj and esincj). The orbital param- 
eters are given in Table [141 ^^'^ the radial velocity dataset 
is plotted in Figure [9] with residuals shown for a circular 



orbit. The Figure also shows models of a circular and an 
eccentric orbit (with e — 0.023), but they are almost undis- 
tinguishable. For the circular orbit, we obtained = 18.00, 
yielding a value of BICc = 160.30 and for the eccentric or- 
bit, we obtained x^ ~ 15.91 and a value of BICe = 164.48. 
We repeated the calculations and set ar = 8.45 m s~^ for 
SOPHIE to obtain a reduced x^ of unity for the best-fit 
circular orbit (using the SOPHIE data alone), while we set 
ar = Om s'^ for the HIRES data since an eccentric orbit 
for that dataset alone yields a reduced x^ of 0.56, indicating 
over-fitting. For a circular orbit, we obtained ~ 15.97, 
resulting in a value of BICc = 159.38 and for an eccentric 
orbit, we obtained x^ = 13.88 and a value of BICe ~ 163.56. 
We therefore find that the circular orbit cannot be excluded 
for TrES-2b. Furthermore, we exclude the possibility that 
e > 0.1. 

WASP-2 (new HARPS data) 

WASP-2b is a 0.85 Mj planet on a 2.15 day orbit around 



a Kl star (V=12), first reported by Collier Cameron et al. 
( |2007[ ). We use 8 new HARPS radial velocity measurements 
and 7 of the original 9 SOPHIE measurements (we drop the 
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Parameter 



HIRES, O'Donovan et al. 1 2006 i 



HIRES, SOPHIE, this work 



Centre-of-mass velocity Vq [m s"-'^] 
Orbital eccentricity e 
Argument of periastron lu [°] 
e cos id 
esinw 

Velocity semi-amplitude K [ms^^] 



(adopted) 
(unconstrained) 



181.3±2.6 



-29.8±2.4 (HIRES), -315.5±5.0 (SOPHIE) 
0.023±0.014, e < 0.051 
(unconstrained) 
0.002±0.009 
-0.022±0.016 
181.1±2.5 



Table 14. System parameters for TrES-2. Left: |0'Donovan et al.| | [2006| . Right: Results from our HARPS radial velocity data. Median 
values for Vi) and K are quoted for the circular orbits, as well as 68.3% confidence limits obtained from the eccentric solution. The 95% 
limit on the eccentricity is also given. 
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Figure 9. Plot showing our new SOPHIE radial velocity data for TrES-2, plotted against time (left) and orbital phase with respect to 
Ttr (right). A circular orbit (solid line) and an orbit with the best-fit eccentricity (dotted line, but almost undistinguishable from the 
circular solution since e = 0.023) are overplotted. The residuals relative to the circular orbit are shown in the bottom panel. 



first measurement, which has an uncertainty about 15 times 
larger than the rest, and the fifth, which shows a 3-a" devia- 



tion at a phase close to the transit) in Collier Cameron et al 



(20071 to work out the orbital parameters of WASP-2b. We 
impose the period P = 2.15222144(39) d and mid-transit 
time Ttr = 2453991.51455(17) BJD as given from photome- 



try in Southworth et al. (20101. We used 15 measurements 



in all, and count the two constraints from photometry as 
two additional datapoints (A'' =17) and used fc = 3 for the 
circular orbit (two Vo, one for each dataset, and the semi- 
amplitude K). 

We estimated the timescale of correlated noise for both 
the HARPS and SOPHIE data to be r = 1.5 d, and we 
estimated ar = 10.4 m s~^ for the SOPHIE data and a,. — 
6.45 m s~^ for the HARPS data to obtain a reduced 
of unity for the circular orbit. We repeated this analysis 
with an eccentric orbit fc = 5 (3 degrees of freedom for the 
circular orbit, and 2 additional degrees of freedom for the 
eccentricity, ecosoj and esinu). The orbital parameters are 
given in Table [Ts] and the radial velocity dataset is plotted 
in Figure [To] with residuals shown for a circular orbit. The 
Figure also shows models of a circular and an eccentric orbit 
(with e — 0.027), but they are almost undistinguishable. For 



the circular orbit, we obtained = 15.60, giving a value of 
BICc ~ 115.08 and for the eccentric orbit, we obtained ~ 
13.88 giving a value of BICe = 119.02. We repeated these 
calculations to obtain a reduced x^ of unity for the eccentric 
orbit and estimated ar = 10.4 m s^^ for the SOPHIE data 
(the SOPHIE dataset did not allow the MCMC to converge 
and yield a reduced x^ of unity with an eccentric orbit) and 
ar = 7.05 m s~^ for the HARPS data. For the circular orbit, 
we obtained x^ ~ 15.16, and a value of BICc ~ 115.47 and 
for the eccentric orbit, we obtained x^ ~ 13.49 and a value 
of BICe = 119.47. We therefore find that the circular orbit 
cannot be excluded for WASP-2. Furthermore, we exclude 
the possibility that e > 0.1. 

Other planets 

HD 189733b and HD209458b are both on orbits that are 



compatible with a circular model: Laughlin et al. ( 2005 \ 



reported the 95% limits on eccentricity for HD 209458b 
(e < 0.042) and we estimate the upper limit for HD189733b 
from Triaud et al. ( 2009 1 assuming a Gaussian probability 



distribution, e < 0.008). In both cases, the eccentricity is 
strongly constrained by the timing of the secondary eclipse. 
No radial velocity data was found for Kepler-5 in the lit- 
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Parameter 



SOPHIE, Collier Cameron et al. ( 2007 1 



SOPHIE and HARPS, this work 



Centre-of-mass velocity Vq [m s^^] 
Orbital eccentricity e 
Argument of periastron ui [°] 



ecosij 
e sinoj 



Velocity semi-amplitude K [ms ^] 



-27863±7 
(adopted) 
(unconstrained) 



155±7 



-27862±7.4 (SOPHIE), -27739.81±4.1 (HARPS), 
0.027±0.023 (<0.072) 
(unconstrained) 
-0.003±0.003 
-0.027±0.027 
156.3±2.1 



Table 15. System parameters for WASP-2. Left: [Collier Cameron et al.| | [2b07[ l. Right: Results from our HARPS radial velocity data. 
Median values for Vb and K are quoted for the circular orbits, as well as 68.3% confidence limits obtained from the eccentric solution 
(see section Analysis) and 95% limit on eccentricity. 
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Figure 10. Plot showing our new HARPS radial velocity data for WASP-2, plotted against time (left) and orbital phase with respect 
to Ttr (right). A circular orbit (solid line) and an orbit with the best-fit eccentricity (dotted line, but almost undistinguishablc from the 
circular solution since e = 0.027) are overplotted. The residuals relative to the circular orbit are shown in the bottom panel. 



erature or online, but we include the results of [Kipping fc| 
Bakos (20111 in this study: Kepler-5b has an eccentricity 



CoRoT-9b is a 0.84 Mj planet on a 95.3 day orbit around 



of e 



O.OSAtoMl, with a 95% upper limit of e < 0.086. 



We therefore classify Kepler-5b as having a circular orbit. 
We also omitted an analysis of the two-planet system HAT- 
P-13, choosing to estimate the 95% limits on the orbital 
eccentricity of HAT-P-13b from the literature (e < 0.022) 
and classify this orbit as circular. 

4.3 Planets on eccentric orbits 



In contrast to Section [4.1| in this Section, we confirm the 
eccentricities of 10 planets. We verify the eccentricities of 
CoRoT-9b, GJ-436b and HAT-P-2b as a test for our pro- 
cedures and we also confirm the eccentricities of HAT-P- 
16b and WASP- 14b, with the former being the planet on 
a short period orbit with the smallest confirmed eccentric- 
ity, and the latter being the planet with the shortest period 
orbit having a confirmed eccentricity. Finally we note the 
confirmed orbital eccentricities of CoRoT-lOb, HAT-P-15b, 
HD17156b, HD80606b and XO-3b. 

CoRoT-9 



a G3 star (V=13.5), first reported by Deeg et al 
who found an eccentricity of e 



(20101, 



0.11 ± 0.0 4 . We used the 
14 HARPS measurements from Deeg et al. (20101, setting 
r — 1.5 d and ar = 3.7 m s~ to obtain a value of re- 
duced of unity for the circular orbit. We imposed the 
prior information from photometry P = 95.2738(14) and 
Ttr = 2454603.3447(1) fromjPeeg eTaT] ([20To|) and obtained 



a value of Xc = 14.05 and xi = 7.90. Using N = 16, = 2 
and fce = 4, we obtain BICc = 106.17 and BICe = 105.57, 
which provides marginal support for an eccentric orbit at 
e = 0.111 ± 0.046, with the 95% limit at e < 0.20. We 
repeated the caculations, setting (Tr = m since this 
results in a reduced x^ of less than unity for the eccentric 
orbit. This time, we obtained a value of Xc ~ 16.65 and 
Xe ~ 9.63. Using A'^ — 16, kc = 2 and ke — 4, we obtain 
BICc = 105.66 and BICe = 104.18, which supports an ec- 
centric orbit at e = 0.111 ± 0.039. 



GJ-436 

GJ-436b is a 0.071 Mj planet on a 2.64 day eccentric orbit 
around a M2.5 star (V=10.7), first reported by jButler et al.[ 
(20041. Deming et al. (20071 detected the secondary eclipse 
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using Spitzer, placing a constraint on the secondary eclipse 
phase </>occ = 0.587 ± 0.005. This translates into ecosuj — 
0.1367 ± 0.0012, which we apply as a Bayesian prior in the 
calculation of our merit function. 

We used the 59 HIRES measurements from IManessI 



et al. (20071, setting r = 1.5 d and a,- = 5.5 m 



to obtain a value of reduced of unity for the circular 
orbit. We imposed the prior information from photome- 



try P = 2.64385(9) fro m [Maness et al.| (|2007| )and Ttr = 
2454280.78149(16) from |Deming et al.| ( |2007[ ) and obtained 
a value of xl = 59.74 and xt = 38.86. Using = 61 (59 
measurements and 2 priors from photometry) and fee = 2 
for the circular orbit, we obtain BICc ~ 371.38. Using 
A'^ = 62 (59 measurements and 3 priors from photometry) 
and fee ~ 3 (Vb, K, esincj) for the eccentric orbit, we ob- 
tain BICe ~ 354.67, which supports an eccentric orbit at 
e = 0.157 ± 0.024, with the 95% limit at e < 0.21. We re- 
peated the caculations, setting ar = 3.95 m s~^ to obtain a 
reduced x^ of unity for the eccentric orbit. This time, we ob- 
tained a value of Xc ~ 88.20 and Xe = 59.39. This time, we 
obtain BICc = 372.67 and BICe = 348.03, which supports 
an eccentric orbit at e = 0.153±0.017, which is in agreement 



with Deming et al. (20071, who reported e = 0.150 ± 0.012. 



HAT-P-16 

HAT-P-16b is a 4.19 Mj planet on a 2.78 day orbit around a 



F8 star (V=10.7), first reported by Buchhave et al. (20101 



The original authors found an eccentricity of e = 0.036 ± 
0.004. We re-analysed the 7 high resolution FIES measure- 
ments, 14 medium resolution FIES measurements and 6 
HIRES measurements, with two priors from photometry on 
the period and mid-transit time. We set r = 1.5 d for all 
instruments and set ar — 115, 185, and 28 m s~^ respec- 
tively for the three instruments to obtain a reduced x^ of 
unity for each individually. We then analysed them together 
using both a circular (x^ = 28.83) and an eccentric orbit 
{x^ = 3.81). Using A'^ = 29, kc = 4 and fee = 6, we obtain 
BICc = 314.74 and BICe = 296.45, which supports an ec- 



centric orbit a t e = 0.034 ± 0. 10. F igure [Tl] (left) shows 

( 2010| , with a circular or- 



the data from Buchhave et al. 



bit overplotted with a solid line and an eccentric orbit with 
the dotted line. The residuals are plotted for the circular 
solution and they show a clear periodic signal. 

We repeated the analysis, this time setting ar = 
(reduced x^ = 0.62, indicating over-fitting), 16 (reduced 
X^ = 33), and 4.7 m s~^ respectively and separately for the 
three datasets (i.e. aiming for a reduced x^ of unity for each 
dataset individually, with an eccentric orbit). We then anal- 
ysed them together using both a circular (x^ = 347.86) and 
an eccentric orbit (x^ = 44.62). Using A'' = 29, fee = 4 and 
fee = 6, we obtain BICe = 541.64 and BICe = 245.14, which 
supports an eccentric orbit at e = 0.034 ± 0.003. We thus 
confirm the eccentricity of HAT-P-16b, which means this 
is the planet with the smallest eccentricity that is reliably 
measured. This is in part helped by the fact that HAT-P-16b 
is a very massive planet, making the radial velocity signal 



for an eccentric orbit very clear. Figure 11 (right) shows the 
data from Buchhave et al. (2010j) again, with an eccentric 



orbit overplotted with the dotted line. 
WASP-14 

WASP-14b is a 7.3 Mj planet on a 2.24 day orbit around a 



F5 star (V=9.8), first reported by Joshi et al 
found an eccentricity of e = 0.091 ± 0.003. |f 



(20091, who 



et al.| 



( {2011 1 confirmed the eccentricity of the orbit and updated 
the precise value to e = 0.088 ± 0.003. This makes WASP- 
14b the planet that is closest to its host star but still has an 
eccentric orbit, taking the place from WASP- 12b. 

CoRoT-10, HAT-P-2, HAT-P-15, HD17156, 
HD80606 and XO-3 

The orbits of the planets CoRoT-lOb (e = 0.110 ± 0.039), 
HAT-P-2b (e = 0.517 ± 0.003), HAT-P-15b (e = 
0.190 ± 0.019), HD17156b (e = 0.677 ± 0.003), HD80606b 
(e = 0.934 ± 0.001) and XO-3b (e = 0.287 ± 0.005) are 
clearly eccentric from existing literature (See for example 
Bonomo et al.||2010| [Loeillet et al.||2008l |Kovacs et al.| 



2010[ [Nutzman et al.||2011| [Hebrard et al.||2010| |200'8r 

respectively) . 



4.4 Planets with orbits that have poorly 
constrained eccentricities 

For 26 of the transiting planets that we attempted to place 
upper limits on their eccentricities, we obtained limits that 
were larger than 0.1. We considered these eccentricities to 
be poorly determined. We discuss the cases of HAT-P-4b, 
WASP-7, XO-2b and Kepler-4b below. 

HAT-P-4 (new SOPHIE data) 

HAT-P-4b is a 0.68 Mj planet on a 3.06 day orbit around 



an F star (V=11.2), first reported by |Kovacs et al.| ( |2007 |. 
We use 13 new SOPHIE radial velocity measurements and 



the 9 HIRES measurements in Kovacs et al. ( 2007 1 to work 



out the orbital parameters of HAT-P-4b. We impose the 
period P — 3.056536(57) d and mid-transit time Ttr = 



24542 4 8.871 6(6) BJD as given from photometry in Kovacs 
et al.| ( |2007| ). We set r = 1.5 d and ar = 3.35 m s"^ for 
SOPHIE and ar = 3.75m s"^ for HIRES, to obtain a re- 
duced x^ of unity for each dataset separately for the best-fit 
circular orbit. We used 22 measurements in all, and count 
the two constraints from photometry as two additional dat- 
apoints (A'=24), and used fc = 3 for the circular orbit (two 
Vb, one for each dataset, and the semi-amplitude K). We 
repeated this analysis with an eccentric orbit fc = 5 (three 
degrees of freedom for the circular orbit, and two additional 
degrees of freedom for the eccentricity, ecoscj and esinoj). 
The orbital parameters are given in Table |16[ and the ra- 
dial velocity dataset is plotted in Figure |12[ with residuals 
shown for a circular orbit. The Figure also shows models 
of a circular and an eccentric orbit (with e — 0.064). For 
the circular orbit, we obtained x^ ~ 22.05, giving a value 
of BICe = 161.96 and for the eccentric orbit, we obtained 
X^ = 16.77 giving a value of BICe = 163.04. We repeated 
these calculations by setting r = 1.5d, ar — 1.81 m s~^ for 
HIRES, and kept ar = 3.35 m s"^ for SOPHIE, since we 
were unable to determine a value of ar that would allow 
the MCMC chain to converge and lead to a x^ of unity for 
an eccentric orbit. This time, we obtained x^ = 25.88 for 
the circular orbit, giving a value of BICc ~ 161.96 and for 
the eccentric orbit, we obtained x^ = 20.05 giving a value 
of BICe = 162.49. We find that the circular orbit cannot 
be excluded for HAT-P-4b, but because the eccentricity is 
e = 0.064 ± 0.028 with an upper limit of e < 0.11, which 
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Figure 11. Plot showing radial velocity data from |Buchhave et al.|p010[ | for HAT-P-16, plotted against orbital phase with respect to 
Ttr- Left: A circular orbit is overplotted with a solid line and an eccentric orbit (e = 0.034) is plotted with a dotted line. The bottom 
panel shows the residuals for a circular orbit: these show a clear periodic signal, indicating the possibility of an eccentric orbit. Right: 
An eccentric orbit (e = 0.034) is plotted with a solid line. The residuals are shown for the eccentric orbit. 



Parameter 



HIRES. 



Kovacs et al. 



1 2007 1 



HIRES+SOPHIE, this work 



Centre-of-mass velocity Vq [m s" 
Orbital eccentricity e 
Argument of pcriastron u [°] 
e cos io 
e sin Lo 

Velocity semi-amplitude K [ms~ 



12.1±0.9 

(adopted) 
(unconstrained) 



81.1±1.9 



20.3 ± 2.6 (HIRES), -1402.0 ±4.0 (SOPHIE) 
0.064±0.028, e < 0.11 
(unconstrained) 
-0.018±0.012 
-0.061±0.027 
81.3±2.6 



Table 16. System parameters for HAT-P-4. Left: p<ovacs et al.| | |2007l l. Right: Results from our new SOPHIE radial velocity data and 
the original HIRES data. Median values for Vb aud K are quoted for the circular orbits, as well as 68.3% confidence limits obtained from 
the eccentric solution. The 95% upper limit on eccentricity is also given. 



is above 0.1, we classify HAT-P-4b as having a poorly con- 
strained eccentricity. 

WASP-7 (new HARPS data) 



a F5 star (V=9.5), first reported by |Hellier eFaL] ( [2009bl ). 
We analysed our 11 new H[ARPS measurements for WASP- 



7 as well as 11 measurements from Hellier et al. (2009b I 



using CORALIE, and used the photometric constraints on 
the orbital period P = 4.954658(55) and mid-transit time 
Ttr = 2453985.0149(12) from the same paper. For both in- 
struments, we set r = 1.5 d and for CORALIE, we set 
(jr = 28.3 m s"^ while for HARPS, we set Or = 210 m s ^ in 
order to get a value of reduced equal to unity for the circu- 
lar orbit. We performed the MCMC analysis twice: the first 
time fitting for the systemic velocity «o and semi-amplitude 
K, and the second time adding two parameters ecosoj and 
esino; to allow for an eccentric orbit. The best-fit parame- 



ters are given in Table 17 1. We plot the radial velocity data 



against time (Figure 13 left) and phase (Figure 13 right). 
When the residuals for a circular orbit are plotted, and a 
scatter of about 30 m s~^ is clearly seen, which is much 



larger than the median uncertainties of cr = 2.21 m s ^ 
on the radial velocity measurements. This is similar to that 



found by [Hellier et al.| (|2009b| from their CORALIE data. 
An eccentric orbit does not reduce the scatter. The value of 

for the circular orbit is 22.37 and that for an eccentric 
orbit is 18.11. This leads to a value of BICc = 250.62 and 
BICe = 252.72, respectively, for 22 measurements, 2 con- 
straints from photometry and 3 and 5 free parameters re- 
spectively (Keplerian orbits, but with two Vb to account for 
a possible offset between the two instruments). This shows 
that the circular orbit is still preferred, and an eccentric 
orbit does not explain the scatter. We repeated this using 
(Jr = 33.8 m s,-^ for CORALIE while for HARPS, we set 
(Jr = 158.5 m s~^ in order to get a value of reduced equal 
to unity for the eccentric orbit. We performed the MCMC 
analysis both for a circular and eccentric orbit. The value of 

for the circular orbit is 146.95 and that for an eccentric 
orbit is 146.86. This leads to a value of BICc = 367.47 and 
BICe = 373.91, respectively, for 22 measurements, 2 con- 
straints from photometry and 3 and 5 free parameters re- 
spectively (Keplerian orbits, but with two Vq to account for 
a possible offset between the two instruments). This shows 
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Figure 12. Plot showing our new HARPS radial velocity data for HAT-P-4, plotted against time (left) and orbital phase with respect to 
Ttr (right). A circular orbit (solid line) and an orbit with the best-fit eccentricity (dotted line, e = 0.064) are overplotted. The residuals 
relative to the circular orbit are shown in the bottom panel. 



that the circular orbit is stiU preferred, and an eccentric or- 
bit does not explain the scatter. WASP-7 is an F5V star, 
with a temperature of TcS = 6400 ± 100 K. Despite the 
result of the original paper that WASP-7 is not chromo- 
spherically active above the 0.02 mag level, [Lagrange et al.| 
1 2009 1 found evidence for other F5V stars showing radial ve- 



locity variability with a scatter at this level, for example HD 
111998, HD 197692 or HD 205289, with scatters of 40 m s"^ 
30 m s~^ and 29 m s~^ respectively. Our derived value of 
eccentricity is e = 0.103 ± 0.061, with the 95% upper limit 
is at e < 0.25. We therefore classify the eccentricity of the 
orbit of WASP-7b as poorly constrained 

In Figure [Ts] we have also plotted the bisector span, the 
signal to noise at order 49, the contrast and full width at half 
maximum for the cross-correlation function against the same 
time axis. The large scatter in radial velocity residuals can 
be seen to be correlated with both the bisector span and the 
full width at half maximum of the cross correlation function. 

XO-2 (new SOPHIE data) 

XO-2 is a 0.6 Mj planet on a 2.62 day orbit around a KO star 
(V=11.2), first reported by|Burke et"ar|(|2007|). We use 9 new 



SOPHIE radial velocity measurements and the 10 HJS mea- 
surements in Burke et al. ( 2007 1 to work out the orbital pa- 
rameters of XO-2. We impose the period P = 2.6158640(21) 
d and mid-transit time Ttr = 2454466.88467(17) BJD as 

We set 



(20091 



given from photometry in Fernandez et al. 
r = 1.5 d and ar = 5.3 m s"^ for SOPHIE and ar = 
m for HJS (because the HJS data alone, with a circu- 
lar orbit, yield a reduced oi 0.78, indicating overfitting) 
to obtain a reduced of unity for the best-fit circular or- 
bit. We used 19 measurements in all, and count the two 
constraints from photometry as two additional datapoints 
(A'^ = 21), and used fc = 3 for the circular orbit (two Vb, one 
for each dataset, and the semi-amplitude K). We repeated 
this analysis with an eccentric orbit k = 5 (two degrees of 



freedom for the circular orbit, and two additional degrees of 
freedom for the eccentricity, e cos lu and e sin cj) . The orbital 
parameters are given in Table |18| and the radial velocity 
dataset is plotted in Figure |14| with residuals shown for a 
circular orbit. The Figure also shows models of a circular 
and an eccentric orbit (with e — 0.064). For the circular or- 
bit, we obtained = 19.65, giving a value of BICc = 165.55 
and for the eccentric orbit, we obtained x^ ~ 17.57 giving 
a value of BICe = 169.55. We repeated the calculations us- 
ing (Tr — 7.05 m s~^ for SOPHIE and ar — ra s~^ for 
HJS (because the HJS data alone, with an eccentric orbit, 
yield a reduced x^ of 0.56, indicating overfitting) to obtain 
a reduced x^ of unity for the best-fit eccentric orbit. For 
the circular orbit, we obtained x^ = 18-02, giving a value 
of BICc = 165.21 and for the eccentric orbit, we obtained 
X^ = 16.01 giving a value of BICe = 169.29. In both cases, 
i.e. using the optimal value of ar for a circular orbit and us- 
ing the optimal value of ar for an eccentric orbit, a circular 
orbit is favoured. The 95% upper limit is e < 0.14, which is 
above 0.1, so we classify the orbital eccentricity of XO-2 as 
poorly constrained. 



0.25+S li, with 



KepIer-4 

Kepler-4b has a derived eccentricity of e = u.zo_q j^ji 

a 95% upper limit of e < 0.43 ( [Kipping fc Bakos||2011[ ), so 
we classify it as "poorly constrained eccentricity". 

Other objects 

For the 8 objects CoRoT-6, HAT-P-1, HAT-P-3, HAT-P- 
6, HD149026, Kepler-6, WASP-10 and WASP-21, we found 
the BICe for an eccentric orbit was smaller than the BICc 
for a circular orbit if we assume a cr,- that yields a reduced 
X^ of unity for an eccentric orbit, whereas the BICc for a 
circular orbit was smaller than the BICe for an eccentric 
orbit if we assume a Ur that yields a reduced x^ of unity for 
a circular orbit. This suggests that the current RV datasets 
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Figure 13. HARPS measurements of WASP-7 plotted against time (left) and phase with respect to the mid-transit time Ttr (right). 
In each case, a solid line is overplotted to represent a circular orbit and the residuals are plotted for this circular orbit. It is clear that 
a signal is present in the residuals (see text). An eccentric orbit with the best-fit value of e = 0.103 is overplotted in both panels with a 
dotted line, but it is almost indistinguishable from the circular solution at this scale. 



Parameter 



Hellier et al. 1 2009b I 



HARPS, this work 



Centre-of-mass velocity Vq [m s"-'^] 
Orbital eccentricity e 
Argument of periastron uj [°] 

e cos ui 
esrau) 

Velocity semi-amplitude K [ms^^] 



-29850.6±1.7 
(adopted) 
(unconstrained) 



97±13 



-29455±103 
0.103 ±0.061 (<0.25) 
(unconstrained) 
0.021±0.068 
0.101±0.074 
96±14 



Table 17. System parameters for WASP-7. Left: |Hellier et al.| ( |2009b[ |. Right: Results from our HARPS radial velocity data. Median 
values for Vb and K are quoted for the circular orbits, as well as 68.3% confidence limits obtained from the eccentric solution (see section 
Analysis) . 



do not constrain the orbit enough for us to detect a finite 
eccentricity. We have aheady discussed the case of WASP- 
10b in Section |4JJ above. 



4.5 Additional planetary systems 

In addition to the 64 planets considered so far, we now in- 
clude 3 additional planets on eccentric orbits, 11 planets 
on orbits where e > 0.1 is excluded at the 95% level and 



two brown dwarves. The additional planets on eccentric or- 
bits are HA T-P-17b ([Howard et al.|2010| e = 0.346 ±0.007), 



HAT-P -21b (|Bakos et al.|2011| e = 0T228±0.016) and HAT 



P-31b ( [Kipping et al.||2011b| e = 0.245 ± 0.005). The addi- 
tional planets on orbits that are consistent with circular are 
CoRoT-18b (e < 0.08 at 3-ff, [Hebrard et al.|20lT | 



HAT-P-20b (e < 0.023, estimated from 
HAT-P-22b (e < 0.031, estimated from 
HAT-P-25b (e < 0.068, estimated from 



Bakos et al. 



Bakos et al. 



Quinn et al 



20111, 



20111, 



20121 



23 



Parameter 



HJS, Burke et al. (2007 1 



HJS, SOPHIE, this work 



Centre-of-mass velocity Vq [m s"-"^] 
Orbital eccentricity e 
Argument of periastron ui [°] 
e cosij 
esino; 

Velocity semi-amplitude K [ms"-*^] 



(adopted) 
(unconstrained) 



85±8 



-1.3 ±6.3 (HJS), 46860.1±4.1 (SOPHIE) 
0.064±0.041 (e < 0.14) 
(unconstrained) 
0.007±0.017 
-0.063±0.047 
98.0±4.0 



Table 18. System parameters for XO-2. Left: [Burke et aT] ( |2007l l. Right: Results from our SOPHIE radial velocity data. Median values 
for Vb and K are quoted for the circular orbits, as well as 68.3% confidence limits obtained from the eccentric solution. The 95% upper 
limit on eccentricity is also given. 
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Figure 14. Plot shov^ing our new SOPHIE radial velocity data for XO-2, plotted against orbital phase with respect to Ttr- A circular 
orbit (solid line) and an orbit with the best-fit eccentricity (dotted line, but almost undistinguishable from the circular solution since 
e = nnnn) are overplotted. The residuals relative to the circular orbit are shown in the bottom panel. 



HAT-P-30b (e < 0.074, estim ated from|Johnson et al.|201l[ ), 
WASP-23b (e < 0.062 at 3-a, [Triaud at al.|2011|), 
WASP-34b (e < 0.058, estim ated from|Smalley et al.|2011[ ), 
WASP-43b (e < 0.04 at 3-cr 
WASP-45b (e < 0.095 
WASP-46b (e < 0.065 



Hellier et al.||2011 l 



Anderson et al. 



Anderson et al. 



20111, 



20111 and 



T Bootis b (e < 0.045, estimated from Butler et al. 2006 1 



The tw o brown dwarves are OGLE-TR-122b ( |Pont et al 



2005b 



2005 



0.205 ± 0.008) and OGLE-TR-123b ( |Pont et al. 
In addition to the above 



= 0). In addition to the above, we also consider the 
case of WASP-38 ( [Barros et al.|[20TT| , which has an eccen- 
tricity of e = 0.031 ± 0.005, indicating it is in the process of 
circularisation, just like WASP-14 and HAT-P-16. 
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5 DISCUSSION 

The Mass-Period plane 

We now discuss the resuhs of the previous sections in the 
context of tidal evolution in hot Jupiters. Figure [15] shows a 
plot of the mass ratio Mp/Ms against orbital period for tran- 
siting planets with orbital period P < 20 days. The empty 
symbols represent orbits that are consistent with circular, 
and the black symbols represent eccentric orbits, whereas 
grey symbols represent objects with small (e < 0.1), but 
significant eccentricities. The circles represent the G dwarfs 
and the squares represent F dwarfs. It appears that the low 
mass hot Jupiters on orbits that are consistent with circu- 
lar around G dwarfs migrate inwards until they stop at a 
minimum period for a given mass, conglomerating on the 



meiss-period relation of Mazeh et al. (2005 1. In this case, the 



heavier planets can move in further before they are stopped. 
Planets heavier than about 1.2 Mj can migrate inwards and 
raise tides on the star, leading to a spin-up of the host star, 
and even synchronisation in some cases where enough angu- 
lar momentum can be transferred from the orbital motion 
into the stellar rotation. In cases where the planetary an- 
gular momentum is insufficient, the process can lead to a 
run-away migration until the planet is destroyed inside the 
star. 

The Roche limit for a planet is defined by Rp — 
0A62aii{Mp/Ms)~'^ . If we write the stopping distance a — 
auR, Ford & Rasio (20061 argued that slow migration on 
quasi-circular orbits would result in a value of a = 1, with 
the only surviving planets being those that stop at their 
Roche limit. On the other hand, if the planets were brought 
in on an eccentric orbit (eg: dynamical interactions within a 
system or capture from interstellar space), and then circu- 
larised by tidal interaction, the value of a should be two. In 
Figure [T5j the dashed line shows this case, with a — 2. This 
does not appear to be a very good fit for the hot Jupiters 
that are on orbits consistent with circular. The dotted lines 
show the range a = 2.5-4.5. As mentionned in |Pont et al.| 
( |2011[ ), this larger value of a could indicate the planets had 
larger radii at the time their orbits were circularised. Sub- 
sequent thermal evolution of the planets would have shrunk 
them (eg: |Baraffe et al.|2004[ ), leaving them further out from 
their current Roche limits. 

Circularisation Timescales 

The process of tidal circularisation, spin-orbit alignment and 
synchronisation are expected to occur roughly in this or- 
der, and over a similar timescale. For close-in systems, this 
timescale is expected to be small compared to the lifetime 
of the system. Hut ( 1981 1 derived equations for the tidal 



evolution due to the equilibrium tide using the assumption 
of weak friction, and constant time-lag At. |Leconte et al.| 
(20101 re-visited this model and showed that the orbital ec- 



centricity evolves according to 



1 de 
e dt 



11 



{Kp [^l.{e)xp^ - if iV,(e)] (11) 
+Ks [ne(e)a;.^-ifiV,(e)]}, 



where Qe{e) and Ne{e) are functions of e and approximately 
equal to unity for small e; Xp and Xs are the cosines of the 



angle between the orbital plane and the planet and stellar 
equators respectively, ojp and ojs are the angular frequencies 
of rotation of the planet and star, and the two terms 



Kp 



and 



Rp 



Ms 
M„ 



f GM^s 



Rs 



/Mpy /R^Y 

\Ms) { a J 



(12) 



(13) 



describe the effect of tides on the planet by the star, and 
vice-versa, respectively, n is the mean orbital motion and 
the semi-major axis is denoted a. Under the assumption of 
a constant-time delay between the exciting tidal potential 
and the response of the equilibrium tide in the relevant body, 
k2,pAtp and fc2,s At^ are constants where k2 are the potential 
Love numbers of degree 2 and At are the constant time lags 
in each of the two bodies. 

We now consider two limits, firstly the case where only 
the tides in the planet dominate, and then the case where 
only tides in the star dominate. When tides in the planet 
dominate, Ks ~ so that we obtain a timescale 



'1 deV 
.e dt/ 



Mr, 



21G fca.pAtp Mi Rl 



(14) 



where we have assumed that — ~ 1, i.e. the 
equation is valid to lowest order in e; uip/n ~ 1, i.e. syn- 
chronisation of the planetary rotation with the orbit and 
Xp ~ 1, i.e. the planet's equator coincides with the orbital 
plane. A similar equation can be written for tides in the 
star, even though Ws/n is not typically unity. As long as 
cjs/n < 18/11, for small e, the effect of tides in the star will 
lead to a decrease in orbital eccentricity. We can therefore 
write, 



Ve dty 



1 



Ms 



21Gk2.sAts 



M2 



Rl 



(15) 



0.01 s and 



We take some typical values of fe2,p^i-p 
fc2,s Ats ~ 1 s, which would co rrespond to tidal quality fac- 
tors ( [Goldreich fc Soter||l966| of about lO'* and lO"* respec- 
tively, in the constant-Q model (in contrast to the constant 
At model that we consider here) for an orbital period of 
about 5 d. 

We expect planets that are further out to be only weak 
affected by tides, whereas close-in planets will experience 
strong tides. Some of these close-in planets will be heavy 
enough and close enough to exert their own infiuence on the 
star by raising stellar tides. This can be seen in Figure \T6\ 
where we have plotted the timescale of circularisation as- 
suming tides inside the star alone against the timescale of 
circularisation assuming tides in the planet alone. The open 
symbols represent orbits that are consistent with circular, 
and the black symbols represent eccentric orbits, whereas 
the grey symbols represent objects with small (e < 0.1), but 
significant eccentricities. The dashed lines represent lines of 
constant circularisation timescale, at 1 Myr, 10 Myr, 100 
Myr, 1 Gyr and 10 Gyr. For the G dwarfs, orbits that are 
consistent with circular and eccentric orbits are cleanly seg- 
regated by the 10 Gyr isochrone, with HAT-P-16b (e — 
0.034 ± 0.003) caught in the process of circularisation. For 
the F dwarfs (open symbols), WASP- 14b {T^s = 6475 ± 100 
K) has a small eccentricity e = 0.008 ± 0.003 and XO-3b 
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Figure 15. Plot showing tlie mass ratio v/s period plane, for transiting planets with orbital period P < 20 days. The low mass hot 
Jupiters on orbits that are consistent with circular around G dwarfs migrate inwards until they stop at a minimum period for a given mass, 
conglomerating on the mass-period relation of |Mazeh et al.| ( |2005[ l. The heavier planets can move in towards the star, and synchronise 
their rotations, as CoRoT-3b and r Bootis b did, or if they lack the angular momentum to synchronise the star, they can continue 
migrating inwards towards their destruction, as WASP-18b appears to be doing. The labelled symbols (except for WASP-18b) represent 
objects on orbits that are consistent with circular where the host star rotation is significantly faster than the expected rotation from the 
isochrones of Strassmeier & Hall (1988). Five objects have been marked with a + symbol to mark objects with upper limits greater than 
e < 0.05 that are described in Section |4.5| The dashed line represents a = 2 for Rp = 1.2Rj, while the dotted lines represent a value in 
the range a = 2.5-4.5 in the equation a = aan (see text). The solid line represents a circularisation isochrone at 1 Gyr for tides in the 
planet alone. 



(Teff = 6429 ± 100 K) has an eccentricity of 0.287 ± 0.005, 
whereas CoRoT-3b {T^a = 6740 ± 140 K) is on an orbit that 
is consistent with circular. This suggests that in the dissipa- 
tion factor in hotter stars may vary in an unknown fashion, 
although the small eccentricity of WASP- 14b and the mod- 
erately small eccentricity of XO-3, together with the short 
timescale for stellar tides indicate that tides in the star are 
clearly important even in these cases. 

Hot Neptunes 

GJ-436b is a planet on an eccentric orbit (e = 0.153±0.017) 
in a region of the mass-scale plane where tidal effects on 
the planet are expected to be significant. The planet is a 
hot Neptune so it is possible that the structure is different 
enough that the tidal quality factor Q is very much higher, 
leading to a longer circularisation timescale. In this case, GJ- 
436b would simply not have had enough time to circularise 
its orbit. Another possibility that was initially suggested by 
[Maness et al.| ( |2007| , is that a second companion may be 
present in the system and is pumping up the eccentricity of 



GJ-436b by secular interactions. Further measurements with 
radial velocity ( Ribas et al.|2009 | and photometry (Ballard 
et al.|2010 l appear to rule this possibility out. 



Synchronisation 

Tidal dissipation leading to orbital circularisation can 
occur in either the planet, the star, or both, according to 
the timescale for each case. On the other hand, synchronisa- 
tion of the host star rotation with the orbital motion would 
depend on tidal effects inside the star alone. This would oc- 
cur on a similar timescale as circularisation in the case of 
dissipation in the star alone. Figure [17] shows the same axes 
as Figure |16[ but on the left panel, the red star symbols 
represent objects with excess stellar rotation. In the case 
of CoRoT-3b and r Bootis b, the rotation of the host star 
has been synchronised with the orbital period. Pont (2009) 
also pointed out that HD 189733 and CoRoT-2b were ro- 
tating faster than expected from the isochrones of Strass- 
meier & Hall (1988), even if the stellar rotations were not 
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Figure 16. Plot showing the timescale of circularisation assuming tides inside the star alone (vertical axis) against the timescale of 
circularisation assuming tides in the planet alone (horizontal axis). The dotted lines represent lines of constant circularisation timescale. 
For the G dwarfs (circles), orbits that are consistent with circular and eccentric orbits are cleanly segregated by the 10 Gyr isochrone, with 
HAT-P-16b (e = 0.034 it 0.003) caught in the process of circularisation. For the F dwarfs (squares), WASP-14b has a small eccentricity 
e = 0.008 it 0.003 and XO-3b has an eccentricity of 0.287 it 0.005, whereas CoRoT-3b is on an orbit that is consistent with circular. The 
short timescale for tides in the star, coupled with the relatively small eccentricities of WASP-14 and XO-3, suggest that tidal effects in 
the star are still operating. 



synchronised. We can now confirm that four more objects 
are clearly in this regime: CoRoT-18, HAT-P-20, WASP- 
19 and WASP-43. The rotation periods of these stars and 
the expected rotation periods are shown in Table [20] From 
Figure |17[ we note that the estimated timescale for orbital 
circularisation due to tidal effects in the star alone is less 
than 5 Gyr for the objects WASP-19, WASP-43, CoRoT-2, 
CoRoT-18 and CoRoT-3. This means that tidal dissipation 
in the star could lead to the excess rotation well within the 
lifetime of these stars. On the other hand, the two objects r 
Bootis b and HAT-P-20 have timescales ~ 10 Gyr, while 



HD 189733b has Ts ~ 80 Gyr. Even in this case, it should 
be noted that the tidal dissipation strength would have to 
be stronger by a single order of magnitude for these objects 
to have been spun up by tidal dissipation inside the star. 
Given that the tidal time lag is uncertain by up to about 
two orders of magnitude, this does not sound implausible. 
In contrast, orbital circularisation in many of these cases 
may well have occured due to dissipation in the planet in- 
stead. Planets that are unable to spin-up their parent stars 
to synchronisation may be doomed to destruction. |Hellier| 
et al.| ( ,2009bj pointed out that the existence of WASP-18 at 
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Name 


Prot (d) 


Expected Prot (d) 


CoRoT-2 


4.52±0.02 


36 


CoRoT-18 


6.3±0.9 


49 


HAT-P-20 


11.3±2.2 


57 


HD 189733 


12.95±0.01 


57 


WASP-19 


10.5±0.2 


42 


WASP-43 


7.6±0.7 


57 



Table 20. Table showing the systems with excess rotation in the 
left panel of Figure |17| P^ot stellar rotation period today, 

and 'Expected Pj-ot' expected rotation period of the star 

as estimated from the rotation isochrones of Strassmeier & Hall 
(1988). 



its current position in the mass-period plane suggests that 
either the tidal dissipation in the system is several orders 
of magnitude smaller than expected, or that the system is 
caught at a very special time while it is in the last 10~* of 
the estimated lifetime of the system. The latter possibility 
sounds more plausible, considering the striking paucity of 
heavy planets at short period. 

Spin-orbit alignment 

The right panel of Figure [17] shows the same axes 
(timescales) , but now the circles represent G stars and the 
squares represent F stars. The empty symbols represent 
aligned systems (A < 30°), and the filled symbols represent 
misaligned systems (A > 30°). In this case, the G dwarfs 
are aligned, except for CoRoT-l(Tcff = 5950 ± 150 K) and 
WASP-1 (Tcff = 6110 ± 245 K) are actually hot stars, and 
WASP-8 is outside the region of strong tides in the star. 
CoRoT-1, WASP-1 and the F dwarfs, display a spread in 
terms of aligned and misaligned, even in cases of strong 



tides. Winn et al. (20101 found a link between the pres- 
ence of a convective core and spin-orbit alignment by tidal 
effects. Thus, exoplanets could migrate inwards by planet- 
planet scattering, giving rise to orbits with a range of ec- 
centricities and spin-orbit angles. Planets in orbit around 
cooler stars (Toff < 6250 K, where the stellar convective 
region is significant), can have their orbital angular momen- 
tum aligned with the stellar rotation, while planets in orbit 
around hot stars (Tcff > 6250 K, where the extent of the 
convective region is negligible) manage to keep their initial 
misalignment. 



interactions in close-in exoplanet systems. The present ob- 
servations support a scenario where low mass hot Jupiters 
migrate inwards and circularise their orbits until they stop 
at a minimum period for a given mass, conglomerating on 



the mass-period relation of Mazeh et al. ( 2005 \ . The heavier 



planets are able to move further inwards before they stop. 
Planets heavier than about 1.2 Mj can raise tides on the star 
as they migrate inwards, leading to a spin-up of the host star 
(Pont, ^2009) , and even spin-orbit synchronisation in some 
cases where enough angular momentum can be transferred 
from the orbital motion into the stellar rotation. This ap- 
pears to be the case for CoRoT-3b, r Bootis b, HD 189733, 
CoRoT-2b, CoRoT-18, HAT-P-20, WASP-19 and WASP- 
43, where the first two are synchronised, and the rest show 
clear evidence of excess rotational angular momentum in the 
star. If the planetary angular momentum is insufficient, the 
process can lead to a run-away migration and the planet is 



destroyed, as appears to be the case for WASP- 18b (Hellier 
et al.|200"9 1. This is also supported by the lack of such heavy 
planets at short period. As suggested by Winn et al. (20101, 



tidal effects in G dwarfs are also responsible for aligning the 
spin of the star with the orbit of the planet, whereas the 
same effect is much less effective in the case of the hotter 
F stars. Overall, therefore, the present data on close-in exo- 
planets support the case for a prominent role for tidal inter- 
actions between the planet and the host star in the orbital 
evolution of hot Jupiters. 
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6 CONCLUSION 

We have recalculated estimates of orbital eccentricity for a 
population of known transiting planets and included a noise 
treatment to account for systematic effects in the data. As 



Laughlin et al. ( 2005 1 showed using synthetic data, analysis 
of radial velocity data can result in a derived eccentricity 
at a few a level even in cases where the orbit is in fact 
consistent with circular. In a similar way, correlated noise in 
the instrument or atmosphere, stellar activity, or additional 
companions to the host star can cause a spurious eccentricity 
detection, the cases of WASP-12 and WASP-10 being two 
examples highlighted in this paper. 

Once these confusing effects are accounted for, a much 
clearer picture emerges, highlighting the importance of tidal 
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Figure 17. The axes are the same as Figure [Tg] Left: The star symbols represent objects where there is evidence of spin-up. These are 
stars that rotate faster than predicted by the isochrones of Strassmeier & Hall (1988). In the case of the two hot stars CoRoT-3 and t 
Bootis b, the stellar rotation have even become synchronised with the orbital period. No objects with a stellar tidal dissipation timescale 
larger than about Ts > 10^^ years show any evidence of excess rotation, supporting the case for tidal involvement in the objects with 
excess rotation. 

Right: The circles indicate aligned systems (A < 30°), whereas the star symbols represent misaligned systems (A > 30°). In this case, 
the G dwarfs are aligned (CoRoT-1 and WASP-1 are actually hot stars, and WASP-8 is outside the region of strong tides in the star). 
The F dwarfs, on the other hand, display a spread in terms of aligned and misaligned, even in cases of strong tides, in agreement with 
Winn et al. (2010). 



30 A^. Husnoo et al. 



REFERENCES 

Alonso R., Alapini A., Aigrain S., Auvergne M., Baglin A., 
Barbieri M., Barge P., Bonomo A. S., et al., 2009, A & A, 
506, 353 

Alonso R., Guillot T., Mazeh T., Aigrain S., Alapini A., 

Barge P., Hatzcs A., Pont P., 2009, A & A, 501, L23 
Anderson D. R., Collier Cameron A., Gillon M., Hellier C, 

Jehin E., Lendl M., Maxted P. F. L., Queloz D., et al., 

2011, MNRAS (submitted), arXiv:1105.3179 
Anderson D. R., Gillon M., Hellier C., Maxted P. F. L., 

Pepe P., et al., 2008, MNRAS, L4, 387 
Anderson D. R., Hellier C, Gillon M., Triaud A. H. M. J., 

Smalley B., Hebb L., Collier Cameron A., Maxted P. F. L., 

ot al., 2010, ApJ, 709, 159 
Bakos G. A., Hartman J., Torres G., Latham D. W., Kovacs 

G., Noyes R. W., Fischer D. A., Johnson J. A., et al., 2011, 

ApJ, 742, 116 

Ballard S., Christiansen J. L., Charbonneau D., Deming 

D., Holman M. J., Fabrycky D., A'Hearn M. F., Wellnitz 

D. D., ct al., 2010, ApJ, 716, 1047 
Baraffe I., Selsis P., Chabrier G., Barman T. S., Allard P., 

Hauschildt P. H., Lammer H., 2004, A & A, 419, L13 
Barros S. C. C, Faedi P., Collier Cameron A., Lister T. A., 

McCormac J., PoUacco D., Simpson E. K., Smalley B., 

et al., 2011, A & A, 525, A54+ 
Bonomo A. S., Santerne A., Alonso R., Gazzano J.-C, 

Havel M., Aigrain S., Auvergne M., Baglin A., et al., 2010, 

A & A, 520, A65 
Buchhave L. A., Bakos G. A., Hartman J. D., Torres G., 

Kovacs G., Latham D. W., Noyes R. W., Esquerdo G. A., 

et al., 2010, ApJ, 720, 1118 
Burke C. J., McCullough P. R., Valenti J. A., Johns-KruU 

C. M., Janes K. A., Heasley J. N., et al., 2007, APJ, 671, 

2115 

Butler R. P., Vogt S. S., Marcy G. W., Fischer D. A., 
Wright J. T., Henry G. W., Laughlin G., Lissauer J. J., 
2004, ApJ, 617, 580 

Butler R. P., Wright J. T., Marcy G. W., Fischer D. A., 
Vogt S. S., Tinney C. G., Jones H. R. A., Carter B. D., 
Johnson J. A., McCarthy C, Penny A. J., 2006, ApJ, 646, 
505 

Campo C. J., Harrington J., Hardy R. A., Stevenson K. B., 

Nymeyer S., Ragozzine D., Lust N. B., Anderson D. R., 

et al., 2011, ApJ, 727, 125 
Christian D. J., Gibson N. P., Simpson E. K., Street R. A., 

Skillen I., PoUacco D., Collier Cameron A., Joshi Y. C, 

et al., 2009, MNRAS, 392, 1585 
Collier Cameron A., Bouchy P., Hcbrard G., Maxted P., 

PoUacco D., Pont F., Skillen I., Smalley B., et al., 2007, 

MNRAS, 375, 951 
Deeg H. J., Moutou C, Erikson A., Csizmadia S., Tingley 

B., Barge P., Bruntt H., Havel M., et al., 2010, Nature, 

464, 384 

Deming D., Harrington J., Laughlin G., Seager S., Navarro 
S. B., Bowman W. C, Horning K., 2007, ApJ, 667, L199 
Fernandez J. M., Holman M. J., Winn J. N., Torres G., 
Shporer A., Mazeh T., Esquerdo G. A., Everett M. E., 
2009, AJ, 137, 4911 
Ford E. B., 2006, ApJ, 642, 505 
Ford E. B., Rasio P. A., 2006, ApJ, 638, L45 
Gillon M., Anderson D. R., Triaud A. H. M. J., Hellier C, 



Maxted P. F. L., PoUaeo D., Queloz D., Smalley B., et al., 

2009, A & A, 501, 785 
Gillon M., Smalley B., Hebb L., Anderson D. R., Triaud A. 

H. M. J., Hellier C, Maxted P. F. L., Queloz D., Wilson 

D. M., 2009, A & A, 496, 259 
Goldreich P., Soter S., 1966, Icarus, 5, 375 
Hebb L., Collier-Cameron A., Loeillet B., Pollacco D., He- 

brard G., Street R. A., Bouchy F., Stempels H. C, et al., 

2009, ApJ, 693, 1920 
Hebb L., Collier-Cameron A., Triaud A. H. M. J., Lister 

T. A., Smalley B., Maxted P. F. L., Hellier C, Anderson 

D. R., et al., 2010, ApJ, 708, 224 

Hebrard G., Bouchy P., Pont P., Loeillet B., Rabus M., 

et al., 2008, A & A, 488, 763 
Hcbrard G., Desert J.-M., Diaz R. P., Boisse I., Bouchy 

P., Lecavelier Des Etangs A., Moutou C, Ehrenreich D., 

et al., 2010, A & A, 516, A95+ 
Hebrard G., Evans T. M., Alonso R., Fridlund M., Ofir A., 

Aigrain S., Guillot T., Almenara J. M., et al., 2011, A & 

A, 533, A130 

Hellier C, Anderson D. R., Cameron A. C, Gillon M., 
Hebb L., Maxted P. F. L., Queloz D., Smalley B., et al., 

2009, Nature, 460, 1098 

Hellier C, Anderson D. R., Collier Cameron A., Gillon M., 
Jehin E., Lendl M., Maxted P. F. L., Pepe F., et al., 2011, 
A & A, 535, L7 

Hellier C, Anderson D. R., Gillon M., Lister T. A., Maxted 
P. F. L., Queloz D., et al., 2009, The Astrophysical Journal 
Letters, 690, L89 

Henry G. W., Winn J. N., 2008, AJ, 135, 68 

Holman M. J., Winn J. N., Latham D. W., et al., 2006, 
ApJ, 652, 1715 

Howard A. W., Bakos G. A., Hartman J., Torres G.. Sh- 
porer A., Mazeh T., Kovacs G., Latham D. W., et al., 

2010, arXiv: 1008.3898 

Husnoo N., Pont P., Hcbrard G., Simpson E., Mazeh T., 
Bouchy P., Moutou C, Arnold L., Boisse I., et al., 2011, 
MNRAS, 413, 2500 

Hut P., 1981, A & A, 99, 126 

Johnson J. A., Winn J. N., Bakos G. A., Hartman J. D., 
Morton T. D., Torres G., Kovacs G., Latham D. W., et al., 

2011, ApJ, 735, 24 

Joshi Y. C, Pollacco D., Cameron A. C, Skillen I., Simpson 

E. , Steele I., Street R. A., Stempels H. C, et al., 2009, 
MNRAS, 392, 1532 

Kipping D., Bakos G., 2011, ApJ, 730, 50 

Kipping D. M., Hartman J., Bakos G. A., Torres G., 

Latham D. W., Bayliss D., Kiss L. L., Sato B., et al., 

2011, AJ, 142, 95 
Kovacs G., Bakos G. A., Hartman J. D., Torres G., Noyes 

R. W., Latham D. W., Howard A. W., Fischer D. A., 

et al., 2010, ApJ, 724, 866 
Kovacs G., Bakos G. A., Torres G., Sozzetti A., Latham 

D. W., Noyes R. W., Butler R. P., et al., 2007, The As- 
trophysical Journal Letters, 670, L41 
Lagrange A., Desort M., Galland F., Udry S., Mayor M., 

2009, A & A, 495, 335 
Lanza A. F., Pagano I., Leto G., Messina S., Aigrain S., 

Alonso R., Auvergne M., Baglin A., et al., 2009, AAP, 

493, 193 

Laughlin G., Marcy G. W., Vogt S. S., Fischer D. A., Butler 
R. P., 2005, The Astrophysical Journal Letters, 629, L121 



31 



Leconte J., Chabrier G., Baraffe I., Levrard B., 2010, A & 

A, 516, A64+ 
Liddle A. R., 2007, MNRAS, 377, L74 
Loeillet B., Shporer A., Bouchy F., Pont F., Mazeh T., 

Beuzit J. L., Boisse I., Bonfils X., da Silva R., et al., 2008, 

AAP, 481, 529 

Maciejewski G., Dimitrov D., Neuhauser R., Tetzlaff N., 
Niedzielski A., Raetz S., Chen W. P., Walter F., et al., 

2011, MNRAS, 411, 1204 

Maciejewski G., Errmann R., Raetz S., Seeliger M., Spale- 
niak I., Neuhauser R., 2011, A & A, 528, A65+ 

Maness H. L., Marcy G. W., Ford E. B., Hauschildt P. H., 
Shreve A. T., Basri G. B., Butler R. P., Vogt S. S., 2007, 
PASP, 119, 90 

Mathieu R. D., Mazeh T., 1988, ApJ, 326, 256 

Mazeh T., 2008, in M.-J. Goupil & J.-P. Zahn ed., EAS 
Publications Series Vol. 29 of EAS Publications Series, 
Observational Evidence for Tidal Interaction in Close Bi- 
nary Systems, pp 1-65 

Mazeh T., Zucker S., Pont F., 2005, MNRAS, 356, 955 

Nutzman P., Gilliland R. L., McCuUough P. R., Charbon- 
neau D., Christensen-Dalsgaard J., Kjeldsen H., Nelan 
E. P., Brown T. M., ct al., 2011, ApJ, 726, 3 

O'Donovan F. T., Charbonneau D., Mandushev G., Dun- 
ham E. W., Latham D. W., Torres G., et al., 2006, ApJ, 
651, L61 

Pal A., Bakos G. A., Torres G., Noyes R. W., Latham 
D. W., Kovacs G., Marcy G. W., Fischer D. A., Butler 
R. P., Sasselov D. D., Sipocz B., Esquerdo G. A., Kovacs 
G., Stefanik R., Lazar J., Papp I., Sari P., 2008, ApJ, 680, 

1450 

Pope F., Mayor M., Galland F., Naef D., Queloz D., Santos 

N. C., Udry S., Burnet M., 2002, A & A, 388, 632 
Pont F., 2009, MNRAS, 396, 1789 

Pont F., Bouchy F., Melo C, Santos N. C., Mayor M., 
Queloz D., Udry S., 2005, A & A, 438, 1123 

Pont F., Hebrard G., Irwin J. M., Bouchy F., Moutou C, 
Ehrcnreich D., Guillot T., Aigrain S., et al., 2009, A & A, 
502, 695 

Pont F., Husnoo N., Mazeh T., Fabrycky D., 2011, MN- 
RAS, 414, 1278 

Pont F., Melo C. H. F., Bouchy F., Udry S., Queloz D., 
Mayor M., Santos N. C., 2005, A & A, 433, L21 

Pont F., Zucker S., Queloz D., 2006, MNRAS, 373, 231 

Quinn S. N., Bakos G. A., Hartman J., Torres G., Kovacs 
G., Latham D. W., Noyes R. W., Fischer D. A., et al., 

2012, ApJ, 745, 80 

Raetz S., Mugrauer M., Schmidt T. O. B., Roell T., Eisen- 
bciss T., Hohle M. M., Koeltzsch A., Vaiiko M., Ginski 
C., Marka C., Moualla M., Tetzlaff N., Seifahrt A., Broeg 
C., Koppenhoefcr J., Raetz M., Neuhauser R., 2009, As- 
tronomische Nachrichten, 330, 459 
Rasio F. A., Ford E. B., 1996, Science, 274, 954 
Rauer H., Queloz D., Csizmadia S., Dcleuil M., Alonso R., 
Aigrain S., Almenara J. M., Auvergne M., et al., 2009, A 
& A, 506, 281 

Ribas I., Font-Ribera A., Beaulieu J.-P., Morales J. C, 
Garcfa-Melendo E., 2009, in lAU Symposium Vol. 253 of 
lAU Symposium, The case for a close-in perturber to GJ 
436 b. pp 149-155 

Sivia D., 2006, Data Analysis: A Bayesian Tutorial. Oxford 
University Press 



Smalley B., Anderson D. R., Collier Cameron A., Hel- 
lier C, Lendl M., Maxted P. F. L., Queloz D., Triaud 
A. H. M. J., et al., 2011, A & A, 526, A130-I- 

Southworth J., Hinse T. C, J0rgensen U. G., Dominik M., 
Ricci D., Burgdorf M. J., Hornstrup A., Wheatley P. J., 
et al., 2009, MNRAS, 396, 1023 

Southworth J., Mancini L., Novati S. C, Dominik M., 
Glitrup M., Hinse T. C, J0rgensen U. G., Mathiasen M., 
et al., 2010, MNRAS, 408, 1680 

Triaud A. H. M. J., Queloz D., Bouchy F., Moutou C, 
Cameron A. C, Claret A., Barge P., Benz W., et al., 2009, 
A & A, 506, 377 

Triaud A. H. M. J., Queloz D., Hellier C, Gillon M., Smal- 
ley B., Hcbb L., Collier Cameron A., Anderson D., et al., 
2011, A & A, 531, A24-h 

Welsh W. F., Orosz J. A., Seager S., Fortney J. J., Jenkins 
J., Rowe J. F., Koch D., Borucki W. J., 2010, ApJ, 713, 
L145 

Wilson D. M., Gillon M., Hellier C, Maxted P. F. L., Pepe 
F., et al., 2008, The Astrophysical Journal Letters, 675, 
L113 

Winn J. N., Fabrycky D., Albrecht S., Johnson J. A., 2010, 

ApJ, 718, L145 
Winn J. N., Holman M. J., Carter J. A., Torres G., Osip 

D. J., Beatty T., 2009, AJ, 137, 3826 
Winn J. N., Johnson J. A., Albrecht S., Howard A. W., 

Marcy G. W., Crossfield I. J., Holman M. J., 2009, ApJ, 

703, L99 

Winn J. N., Noyes R. W., Holman M. J., Charbonneau D., 
Ohta Y., Taruya A., Suto Y., Narita N., et al., 2005, ApJ, 
631, 1215 



